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Abstract. Many procedures in science, engineering and medicine produce data 
in the form of geometric shapes. MathematicaUy, a shape can be modeled as 
an un-parameterized immersed sub-manifold, which is the notion of shape used 
here. Endowing shape space with a Riemannian metric opens up the world of 
Riemannian differential geometry with geodesies, gradient flows and curvature. 
Unfortunately, the simplest such metric induces vanishing path-length distance 
on shape space. This discovery by Michor and Mumford was the starting point 
to a quest for stronger, meaningful metrics that should be able to distinguish 
salient features of the shapes. Sobolev metrics are a very promising approach 
to that. They come in two flavors: Outer metrics which arc induced from 
metrics on the diffeomorphism group of ambient space, and inner metrics which 
are defined intrinsically to the shape. In this work, Sobolev inner metrics are 
developed and treated in a very general setting. There are no restrictions on the 
dimension of the immersed space or of the ambient space, and ambient space 
is not required to be flat. It is shown that the path- length distance induced by 
Sobolev inner metrics docs not vanish. The geodesic equation and the conserved 
quantities arising from the symmetries are calculated, and well-posedness of the 
geodesic equation is proven. Finally examples of numerical solutions to the 
geodesic equation are presented. 
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Chapter 1 

Introduction 



1.1 The Riemannian setting for shape analysis 

From very early on, shapes spaces have been analyzed in a Riemannian setting. 
This is also the setting that has been adopted in this work. The Riemannian 
setting is well-suited to shape analysis for several reasons. 

• It formalizes an intuitive notion of similarity of shapes: Shapes that differ 
only by a small deformation are similar to each other. To compare shapes, 
it is thus necessary to measure deformations. This is exactly what is 
accomplished by a Riemannian metric. A Riemannian metric measures 
continuous deformations of shapes, that is, paths in shape space. 

• Riemannian metrics on shape space have been used successfully in com- 
puter vision for a long time, often without any mention of the underlying 
metric. Gradient flows for shape smoothing are an example. An under- 
lying metric is needed for the definition of a gradient. Often, the metric 
that has been used implicitly was the L^-metric that has however turned 
out to be too weak. 

• The exponential map that is induced by a Riemannian metric permits to 
linearize shape space: When shapes are represented as initial velocities of 
geodesies connecting them to a fixed reference shape, one effectively works 
in the linear tangent space over the reference shape. Curvature will play 
an essential role in quantifying the deviation of curved shape space from 
its linearized approximation. 

• The linearization of shape space by the exponential map allows to do 
statistics on shape space. 

A disadvantage of the Riemannian approach is that shapes can be compared 
with each other only when there is a deformation between them. 
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1.2 Related work 

In mathematics and computer vision, shapes have been represented in many 
ways. Point clouds, meshes, level-sets, graphs of a function, currents and mea- 
sures are but some of the possibilities. Furthermore, the resulting shape spaces 
have been endowed with many different metrics. Approaches found in the liter- 
ature include: 

• Inner metrics on shape space of unparametrized immersions. These met- 
rics are induced from metrics on parametrized immersions. Since this is 
the approach studied in this work, more references will be given later. 

• Outer metrics on various shape spaces (images, embedded surfaces, land- 
marks, measures and currents) that the diffeomorphism group of ambient 
space is acting on. See for example [51[T11 . [5g i Hg i [ro i [55 1 [2 ^ [M|. 

• Metamorphosis metrics. See for example [47l [25]. 

• The Wasserstein metric or Monge-Kantorovic metric on shape space of 
probability measures. See for example pi [51 [T^ [T^ . 

• The Weil-Peterson metric on shape space of planar curves. See for example 

021113130]. 

• Current metrics. See for example [48l [19l [20] . 

More references can be found in the review papers [TJ |4j jM] jSTJ |4T1 149] . 

1.3 Inner versus outer metrics 

Outer metrics measure how much ambient space has to be deformed in order to 
yield the desired deformation of the shape. This concept has been introduced by 
the Scottish biologist, mathematician and classics scholar d'Arcy Thompson in 
1942, already. As Thompson declared in the epilogue of his book "On Growth 
and Form" 45 , his aim was to show that "a certain mathematical aspect of 
morphology is essential to (the) proper study of Growth and Form" . In the 
chapter "On the comparison of related forms" of this book, Thompson pictures 
transformations of the ambient space by a Cartesian grid. A transformation of 
the ambient space and the grid then results in a deformation of the embedded 
shape, too. An example is given in figure [TT] 

Thompson's notion of shape transformation and the concept of outer metrics 
is fundamentally different from the notion of shape transformation that underlies 
this work. In this work, so called inner metrics are treated. Inner metrics 
measure how much the shape itself is deformed. A deformation of the shape 
does not carry with it a deformation of the ambient space. The ambient space 
is fixed, and the shape moves within it. This is illustrated in figure [T72j 

Riemannian metrics on shape space measure infinitesimal deformations. An 
infinitesimal deformation of ambient space is a vector field on ambient space and 
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Figure 1.1: Transformation of the ambient space carrying with it a transforma- 
tion from human skull to chimpanzee and baboon. From |45L p. 318-319]. 
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Figure 1.2: The same transformation as in figure |1.1[ but only the shape is 
transformed, not the ambient space. Modification of [lH p. 318-319]. 
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could be pictured as a small arrow attached to every point in ambient space. 
This is the kind of deformation that is measured by outer metrics. In contrast to 
this, an infinitesimal deformation of the shape itself is a (normal or horizontal) 
vector field along the shape. It could be pictured as a small arrow attached to 
every point of the shape. Figure |1.3| illustrates the two kinds of deformation. 




HumtB ikull. Huinta ikuU- 



Figure 1.3: Infinitesimal transformation of ambient space (left) as measured by 
an outer metric and infinitesimal transformation of the shape itself (right) as 
measured by an inner metric. 

The deformations that are measured by inner metrics have much smaller 
(lower dimensional) support than those measured by outer metrics. This is 
an advantage in numerics. A disadvantage is that the differential operator 
governing an inner metric usually depends on the shape, which is not the case 
for outer metrics. 

Section |3.16| explains how the difference between inner and outer metrics 
becomes manifest in the mathematical theory. Sections |3.3| and |3.13| explain in 
mathematical terms what an infinitesimal deformation of a shape is. 



1.4 Where this work comes from 

This work continues the study of inner metrics on shape space of unparametrized 
immersions of a fixed manifold M into a Riemannian manifold N. M fixes the 
topology of the shapes. N is the ambient space, and it is endowed with a fixed 
metric. For example, M could be a simple geometric object like a sphere, and 
N could be some M". 

It came as a big surprise when Michor and Mumford found out in [SHIISZ] that 
the simplest and most natural Riemannian metric on this shape space (the L^- 
metric) induces vanishing geodesic distance on shape space. More precisely the 
result is that any two shapes can be deformed into each other by a deformation 
that is arbitrarily small when measured with respect to this metric. The result 
also holds for L^-metrics on diffeomorphism groups and for outer L^-metrics on 
shape space. The discovery of this degeneracy was the starting point to a quest 
for stronger metrics. These metrics should be able to distinguish salient features 
of the shapes. The meaning of salient of course depends on the application. 

One approach to strengthen the inner metric is by weighting it by a 
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function depending on the mean curvature and/or the volume [551 [571 [711^ [TU] . 
Such metrics have been called almost local metrics. 

Another approach, and the approach taken in this work, is to incorporate 
derivatives of the deformation vector fields in the definition of the metric. This 
yields Sobolev inner metrics on shape space. Sobolev inner metrics have first 
been defined on shape space of planar curves [311 [511 I^H HEl 121] ■ In this work 
and in [5] , these metrics are generalized to higher dimensional shape spaces and 
to a possibly curved ambient space. 

Some parts of this work have been written in collaboration with Martin 
Bauer and can also be found in his thesis IT0\ . This concerns mainly sections [2] 
[3| and |4| covering a lot of background material. 



1.5 Content of this work 

This work progresses from a very general setting to a specific one in three steps. 
In the beginning, a framework for general inner metrics is developed. Then the 
general concepts carry over to more and more specific inner metrics. 

• First, shape space is endowed with a general inner metric, i.e with a 
metric that is induced from a metric on the space of immersions, but 
that is unspecified otherwise. It is shown how various versions of the 
geodesic equation can be expressed using gradients of the metric with 
respect to itself and how conserved quantities arise from symmetries. (This 
is section |3]) 

• Then it is assumed that the inner metric is defined via an elliptic pseudo- 
differential operator. Such a metric will be called a Sobolev-type metric. 
The geodesic equation is formulated in terms of the operator, and existence 
of horizontal paths of immersions within each equivalence class of paths 
is proven. (This is section [s]) Then estimates on the path-length distance 
are derived. Most importantly it is shown that when the operator involves 
high enough powers of the Laplacian, then the metric does not have the 
degeneracy of the L^-metric. (This is section [6]) 

• Motivated by the previous results it is assumed that the elliptic pseudo- 
differential operator is given by the Laplacian and powers of it. Again, the 
geodesic equation is derived. The formulas that are obtained are ready to 
be implemented numerically. (This is section [7j) 

The remaining sections cover the following material: 

• Section [2] treats some differential geometry of surfaces that is needed 
throughout this work. It is also a good reference for the notation that 
is used. The biggest emphasis is on a rigorous treatment of the covariant 
derivative. Some material like the adjoint covariant derivative is not found 
in standard text books. 
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Section |4] contains formulas for the variation of the metric, volume form, 
covariant derivative and Laplacian with respect to the immersion inducing 
them. These formulas are used extensively later. 

Section |8] covers the special case of flat ambient space. The geodesic equa- 
tion is simplified and conserved momenta for the Euclidean motion group 
are calculated. Sobolev-type metrics are compared to the Frechet metric 
which is available in flat ambient space. 

Section [9] treats diffeomorphism groups of compact manifolds as a special 
case of the theory that has been developed so far. 

In section [To] it is shown in some examples that the geodesic equation on 
shape space can be solved numerically. 

Contributions of this work. 

This work is the first to treat Sobolev inner metrics on spaces of immersed 
surfaces and on higher dimensional shape spaces. 

It contains the first description of how the geodesic equation can be for- 
mulated in terms of gradients of the metric with respect to itself when the 
ambient space is not flat. To achieve this, a covariant derivative on some 
bundles over immersions is defined. This covariant derivative is induced 
from the Levi-Civita covariant derivative on ambient space. 

The geodesic equation is formulated in terms of this covariant derivative. 
Well-posedness of the geodesic equation is shown under some regularity 
assumptions that are verified for Sobolev metrics. Well-posedness also 
follows for the geodesic equation on diffeomorphism groups, where this 
result has not yet been obtained in that full generality. 

To derive the geodesic equation, a variational formula for the Laplacian 
operator is developed. The variation is taken with respect to the metric on 
the manifold where the Laplacian is defined. This metric in turn depends 
on the immersion inducing it. 

It is shown that Sobolev inner metrics separate points in shape space when 
the order of the differential operator governing the metric is high enough. 
(The metric needs to be as least as strong as the i/^-metric.) Thus Sobolev 
inner metrics overcome the degeneracy of the L^-metric. 

The path-length distance of Sobolev inner metrics is compared to the 
Frechet distance. It would be desirable to bound Fechet distance by some 
Sobolev distance. This however remains an open problem. 

Finally it is demonstrated in some examples that the geodesic equation for 
the iJ^-metric on shape space of surfaces in can be solved numerically. 
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Differential geometry of 
surfaces and notation 

In this section the differential geometric tools that are needed to deal with 
immersed surfaces are presented and developed. The most important point is a 
rigorous treatment of the covariant derivative and related concepts. 

The notation of [3S] is used. Some of the definitions can also be found in 
[27] . A similar exposition in the same notation is [3 H] . This section has been 
written in collaboration with Martin Bauer and is the same as section 1.1 of his 
Ph.D. thesis JOJ, up to slight modifications. 



2.1 Basic assumptions and conventions 

Assumption. It is always assumed that M and N are connected manifolds of 
finite dimensions m and n, respectively. Furthermore it is assumed that M is 
compact, and that N is endowed with a Riemannian metric g. 



In this work, immersions of M into N will be treated, i.e. smooth func- 
tions M — > N with injective tangent mapping at every point. The set of all 
such immersions will be denoted by Imm(M, N). It is clear that only the case 
dim(M) < dim(A^) is of interest since otherwise Imm(M, N) would be empty. 

Immersions or paths of immersions are usually denoted by /. Vector fields 
on Imm(M, N) or vector fields along / will be called /i, fc, to, for example. Sub- 
scripts like ft = dtf = df /dt denote differentiation with respect to the indicated 
variable, but subscripts are also used to indicate the foot point of a tensor field. 
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2.2 Tensor bundles and tensor fields 

The tensor bundles 

TIM TIM (g) f*TN 



M M 
will be used. Here TJM denotes the bundle of ( )-tensors on M, i.e. 

r s 

TIM ^<^TM®<^ T*M, 

and f*TN is the pullback of the bundle TN via /, see [Ml section 17.5]. A 
tensor field is a section of a tensor bundle. Generally, when E is a bundle, the 
space of its sections will be denoted by T{E). 

To clarify the notation that will be used later, some examples of tensor bun- 
dles and tensor fields are given now. S^T*M = L^y^{TM;R) and A.^T*M = 
L^i^{TM;R) are the bundles of symmetric and alternating (")-tensors, respec- 
tively. n''{M) = r(A'^T*A/) is the space of differential forms, X(M) = T{TM) 
is the space of vector fields, and 

T{f*TN) ^{he C°°{M,TN) : TTjv o /i = /} 

is the space of vector fields along f. 



2.3 Metric on tensor spaces 

Let g € T{S?,qT* N) denote a fixed Riemannian metric on N. The metric 
induced on M by f (z Imm(M, N) is the pullback metric 

g^rge T{SloT*M), g{X, Y) = (/*g)(A, Y) = g{Tf.X, Tf.Y), 

where X, Y are vector fields on M . The dependence of g on the immersion / 
should be kept in mind. Let 

\) = g: TM T*M and (t = 9^^ ■ T*M TM. 

g can be extended to the cotangent bundle T*M = T^M by setting 

5-i(a,/?)=g?(a,/3)=a(/3«) 

for a,/3 G T*M. The product metric 

r s 

extends g to all tensor spaces T^M, and g1®g yields a metric on T^ M (E) f*TN . 
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2.4 Traces 

The trace contracts pairs of vectors and co-vectors in a tensor product: 

Tr : T*M(E) TM = L{TM, TM) ^ Af x M 

A special case of this is the operator ix inserting a vector X into a co- vector or 
into a covariant factor of a tensor product. The inverse of the metric g can be 
used to define a trace 

Tr^ : T*M ® T* M ^ M x M 

contracting pairs of co-vecors. Note that Tr^ depends on the metric whereas Tr 
does not. The fohowing lemma will be useful in many calculations: 

Lemma. 

gl{B, C) ^ Ti{g-^Bg-^C) for B,C e T^M if B or C is symmetric. 
(In the expression under the trace, B and C are seen maps TM — >■ T* M .) 

Proof. Express everything in a local coordinate system u^, . . . u™ of M. 
g^{B, C) = 5^ ( ^ B.kdu' c^^t^ ^ Cjidu^ (g) du'^ 

ik jl 

= Y,g^^B,,g^'C,i = ^g-'^B.^/'Cy = Tv{g-'Bg-'C) 

ijkl ijkl 

Note that only the symmetry of C has been used. □ 

2.5 Volume density 

Let Vol(M) be the density bundle over Af, see [36l section 10.2]. The volume 
density on M induced by / G Imm(A/, N) is 

vol(.g)=vol(r5)er(Vol(A/)). 

The volume of the immersion is given by 

Vol(/) = / vol(/*5) = / vol(5). 

JM Jm 

The integral is well-defined since M is compact. If M is oriented the volume 
density may be identified with a differential form. 

2.6 Metric on tensor fields 

A metric on a space of tensor fields is defined by integrating the appropriate 
metric on the tensor space with respect to the volume density: 

gl{B,C)^ [ gl{Bix),Cix))vo\ig)ix) 
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for B,C e r(T;M), and 

^g{B,C)= f g:®g{B{x),C{x))vQ\{g){x) 

JM 

for S,C e T{TIM ® f*TN), f e Inim(M,7V). The integrals are well-defined 
because M is compact. 

2.7 Covariant derivative 

Covariant derivatives on vector bundles as explained in [311 sections 19.12, 22.9] 
will be used. Let V^, be the Levi-Civita covariant derivatives on {M,g) and 
{N,g), respectively. For any manifold Q and vector field X on Q, one has 

: C°°(Q,rAf) ^ C°°(Q,TM), h^V^h 

W\:C°°{Q,TN)^C°°{Q,TN), h^V^^h. 

Usually the symbol V will be used for all covariant derivatives. It should be 
kept in mind that depends on the metric g = f*g and therefore also on the 
immersion /. The following properties hold [36. section 22.9]: 

1. Vx respects base points, i.e. tt o Vx^ = n o h, where tt is the projection 
of the tangent space onto the base manifold. 

2. Wxh is C°°-linear in X. So for a tangent vector X^ € T^Q, ^ x^h makes 
sense and equals xh){x). 

3. V xh is M- linear in h. 

4. Wx{o,-h) = da{X).h + a.Wxh for a G C°°{Q), the derivation property of 

5. For any manifold Q and smooth mapping q : Q Q and Yy € TyQ one 
has ^Tq.Yyh — ^Yy{h o q)- If F £ X(Qi) and X e X((3) are g-related, 
then S/yih ° q) — iy xh) o q. 

The two covariant derivatives and can be combined to yield a covari- 
ant derivative V x acting on C°°{Q,T^M (g) TN) by additionally requiring the 
following properties section 22.12]: 

6. Vx respects the spaces C°°{Q, T^M ® TN). 

7. V x{h (E) k) — {V xh) (E) k + h (E) (Wxk), a derivation with respect to the 
tensor product. 

8. Vx commutes with any kind of contraction (see [321 section 8.18]). A 
special case of this is 

Vx{a{Y)) = {Vxa){Y) +a{VxY) for a (g) F : iV ^ T^M. 

Property ([I]) is important because it implies that Vx respects spaces of sections 
of bundles. For example, for Q = M and / S C°°(M, N), one gets 

Vx : r(r;M (g f*TN) T{T^M (g) f*TN). 
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2.8 Swapping covariant derivatives 

Some formulas allowing to swap covariant derivatives will be used repeatedly. 
Let / be an immersion, h a vector field along / and X, Y vector fields on M. 
Since V is torsion-free, one has section 22.10]: 

(1) VxTf.Y - VyT/.X - Tf.[X, Y] = Tor(r/.X, Tf.Y) = 0. 
Furthermore one has [3H section 24.5]: 

(2) VxVy/i - Vy^xh - V[x,Y]h = R^o [Tf.X, Tf.Y)h, 
where R^ e fl'^ {N; L{TN,TN)) is the Riemann curvature tensor of {N,g). 

These formulas also hold when / : M x A/ — > iV is a path of immersions, 
h -.M. X M — > TN is a vector field along / and the vector fields are vector fields 
on K X M. A case of special importance is when one of the vector fields is 
(9t,0M) and the other (OR,y), where y is a vector field on M. Since the Lie 
bracket of these vector fields vanishes, ([T]) and ^ yield 

(3) V(a,,OM)r/.(OR,r) - V(o«,y)T/.(9t,0M) = 
and 

(4) V(a,,o„)V(o„y)/i- V(o«,y)V(a„o,,)/^ = i?^(T/.(at, Om), r/.(OK, F))/.. 

2.9 Second and higher covariant derivatives 

When the covariant derivative is seen as a mapping 

V : r(T;M) ^ r(T;+iM) or V : r(r;M(g)/*TA^) ^ r(T;+iM® /*Tiv), 

then the second covariant derivative is simply VV = V^. Since the covariant 
derivative commutes with contractions, can be expressed as 

V| y := iy txV^ = ty Va'V = VxVy - Vy^y for X,Y e X{M). 

Higher covariant derivates are defined according ly as V'^, fc > 0. 

2.10 Adjoint of the covariant derivative 

The covariant derivative 

V : r(r;A/) ^ r(r;+iM) 

admits an adjoint 

V* : r(T;+iM) ^ r(r;M) 

with respect to the metric g, i.e.: 

c^,{^B,C)=gl{B,\/*C). 
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In the same way, V* can be defined when V is acting on r(rjM f*TN). In 
either case it is given by 

V*B = -Trf(VS), 

where the trace is contracting the first two tensor slots of VS. This formula 
will be proven now: 

Proof. The result holds for decomposable tensor fields jS^B £ r(TJ_^j^Af) since 
= [ C^>g:{B,C)vo\{g)~ [ g^V ^,B,C)vo\ig) 

JM JM 

= I -gliB,C)C^>yo\ig)- [ Tr«(/3 ® Vi?), C)vol(g) 

JM JM 

= - div(/3*)B - Tr^(/3 (g) VB), c) 

= - div{/3*)B + Tr5((V/?) <S)B)~ TrS(V(;3 ® S)), c) 

- div(/3»)B + Tr9(V/3)B - Tr3(V(;9 (g) B)), 

= gl(0-TrS(V(/3®i?)),c) 

Here it has been used that Vx5 = 0, that Vx commutes with any kind of 
contraction and acts as a derivation on tensor products |36l section 22.12] and 
that div(X) = Tr(VX) for all vector fields X [36l section 25.12]. To prove the 
resuh for e T{Tl^^M ® f*TN) one simply has to replace gl by gl(^g. □ 

2.11 Laplacian 

The definition of the Laplacian used in this work is the Bochner- Laplacian. It 
can act on all tensor fields B and is defined as 

AS = V*V-B = - Tr^y'^B). 

2.12 Normal bundle 

The normal bundle Nor(/) of an immersion / is a sub-bundle of f*TN whose 
fibers consist of all vectors that are orthogonal to the image of /: 

Nor(/), = {r e r/(,)7V : yX e T,M : g{Y,Tf.X) - O}. 

If dim(M) = dim(iV) then the fibers of the normal bundle are but the zero 
vector. Any vector field h along / € Imm can be decomposed uniquely into 
parts tangential and normal to / as 

h = Tf.h^ + h^, 

where hJ is a vector field on M and is a section of the normal bundle Nor(/). 



2. 1 3. SECOND FUNDAMENTAL FORM AND WEINGARTEN MAPPING13 



2.13 Second fundamental form and Weingarten 
mapping 

Let X and Y be vector fields on M. Then the covariant derivative V xTf.Y 
splits into tangential and a normal parts as 

WxTf.Y = Tf.(VxTf.Y)'^ + (VxT/.r)^ = Tf.VxY + S{X, Y). 

S is the second fundamental form of f. It is a symmetric bilinear form with 
values in the normal bundle of /. When Tf is seen as a section of T*M^ f*TN 
one has S = VTf since 

S{X, Y) = VxTf.Y - Tf.VxY = (VT/)(X, Y). 

The trace of S is the vector valued mean curvature Tr^{S) e r(Nor(/)). 



CHAPTER 2. NOTATION 



Chapter 3 

Shape space 



Briefly said, in this work the word shape means an unparametrized surface. 
(The term surface is used regardless of whether it has dimension two or not.) 
This section is about the infinite dimensional space of all shapes. First an 
overview of the differential calculus that is used is presented. Then some spaces 
of parametrized and unparametrized surfaces are described, and it is shown how 
to define Riemannian metrics on them. The geodesic equation and conserved 
quantities arising from symmetries are derived. 

The agenda that is set out in this section will be pursued in section |5] when 
the arbitrary metric is replaced by a Sobolev-type metric involving a pseudo- 
differential operator and later in section [7] when the pseudo-differential operator 
is replaced by an operator involving powers of the Laplacian. 

This section is common work with Martin Bauer and can also be found in 
section 1.2 of his Ph.D. thesis [lOj . 



3.1 Convenient calculus 

The differential calculus used in this work is convenient calculus [IS]- The 
overview of convenient calculus presented here is taken from [351 Appendix A]. 
Convenient calculus is a generalization of differential calculus to spaces beyond 
Banach and Frechet spaces. Its most useful property for this work is that the 
exponential law holds without any restriction: 

C°° {E xF,G)^ C°° {E, C°° (F, G)) 

for convenient vector spaces E, F, G and a natural convenient vector space struc- 
ture on {F,G). As a consequence variational calculus simply works: For 
example, a smooth curve in C°° {AI, N) can equivalently be seen as a smooth 
mapping M x M. ^ N. The main difficulty in the convenient setting is that 
the composition of linear mappings stops being jointly continuous at the level 
of Banach spaces for any compatible topology. 



15 
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Let E he a locally convex vector space. A curve c : K — > -E is called smooth 
or C°° if all derivatives exist and are continuous - this is a concept without 
problems. Let C°°(]R,£^) be the space of smooth functions. It can be shown 
that C°°(M, E) docs not depend on the locally convex topology of E, but only 
on its associated homology (system of bounded sets). 

E is said to be a convenient vector space if one of the following equivalent 
conditions is satisfied (called c°°-completeness): 

1. For any c e C°°{R,E) the (Riemann-) integral c{t)dt exists in E. 

2. A curve c -.M. ^ E is smooth if and only if A o c is smooth for all A G E' , 
where E' is the dual consisting of all continuous linear functionals on E. 

3. Any Mackey-Cauchy-sequence (i. e. tnm{Xn — Xm) — >■ for some tnm — ^ oo 
in M) converges in E. This is visibly a weak completeness requirement. 

The final topology with respect to all smooth curves is called the c°° -topology 
on E, which then is denoted by c'^E. For Frechet spaces it coincides with the 
given locally convex topology, but on the space V of test functions with compact 
support on M it is strictly finer. 

Let E and F be locally convex vector spaces, and let U C E be c°°-open. 
A mapping f : U ^ F is called smooth or C°°, if / o c e C°"(M,F) for all 
c € C°°(M, [/). The notion of smooth mappings carries over to mappings be- 
tween convenient manifolds, which are manifolds modeled on c°°-open subsets 
of convenient vector spaces. 

Theorem. The main properties of smooth calculus are the following. 

1. For mappings on Frechet spaces this notion of smoothness coincides with 
all other reasonable definitions. Even on this is non-trivial. 

2. Multilinear mappings are smooth if and only if they are bounded. 

3. If f : E ^ U ^ F is smooth then the derivative df : U x E ^ F is smooth, 
and also df : U ^ L{E,F) is smooth where L{E,F) denotes the space of 
all bounded linear mappings with the topology of uniform convergence on 

bounded subsets. 

4-. The chain rule holds. 

5. The space C^{U, F) is again a convenient vector space where the structure 

is given by the obvious injection 



C°°{U,F) 



C°°{R,F)-^ 



Yl C°°(M,R). 



cec°^{R,u),\eF' 



6. The exponential law holds: 



C°°{U, C°°{V, G)) ^C°°{Ux V, G) 



is a linear diffeomorphism of convenient vector spaces. Note that this is 
the main assumption of variational calculus. 



3.2. MANIFOLDS OF IMMERSIONS AND EMBEDDINGS 
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7. A linear mapping f : E ^ C°°{V,G) is smooth (bounded) if and only if 

E — ^ C°°{V,G) G is smooth for each v £ V. This is called the 
smooth uniform boundedness theorem and it is quite applicable. 

Proofs of these statements can be found in l29l . 



3.2 Manifolds of immersions and embeddings 



What has sloppily been called a parametrized surface will now be turned into 
a rigorous definition. Mathematically, parametrized surfaces will be modeled 
as immersions or embeddings of one manifold into another. Immersions and 
embeddings are called parametrized since a change in their parametrization 
(i.e. applying a diffeomorphism on the domain of the function) results in a 
different object. The following sets of functions will be important: 

(1) Emh{M, N) C Imm/(M, N) C Imm(Af, N) C C°°(M, N). 

C°°{M, N) is the set of smooth functions from M to N. Imm(M, N) is the set 
of all immersions of M into N, i.e. all functions / € C°°{M, N) such that T^f 
is injective for all x € M. Imm^(M, TV) is the set of all free immersions. An 
immersion / is called free if the diffeomorphism group of M acts freely on it, 
i.e. f oip ^ f implies ip = Id a/ for all ip £ Diff (M). Emb(M, N) is the set of all 
embeddings of AI into N, i.e. all immersions / that are a homeomorphism onto 
their image. 

The following lemma from 17, 1.3 and 1.4] gives sufficient conditions for an 
immersion to be free. In particular it implies that every embedding is free. 

Lemma. If ip (z Diff(M) has a fixed point and if f o ^ = f for some immersion 
f , then if = Id^v/. 

If for an immersion f there is a point x S f{M) with only one preimage 
then f is free. 



Since M is compact by assumption (see section 2.1 ) it follows that C°°{M, N) 
is a Frechet manifold f55^, section 42.3]. All inclusions in ([T]) are inclusions of 
open subsets: Imm(Af, A^) is open in C°°{M, N) since the condition that the 
differential is injective at every point is an open condition on the one-jet of 
/ [Ml section 5.1]. Imm/(M, A^) is open in Imm(M, A^) by [17, theorem 1.5]. 
Emb(M, A^) is open in Imm/(M, A^) by [29, theorem 44.1]. Therefore all function 
spaces in ([T]) are Frechet manifolds as well. 

When it is clear that M and A^ are the domain and target of the mappings, 
the abbreviations Emb, Imnif, Imm will be used. In most cases, immersions will 
be used since this is the most general setting. Working with free immersions 
instead of immersions makes a difference in section |3.11[ and working with 
embeddings instead of immersions makes a difference in section [6T| The tangent 
and cotangent space to Imm are treated in the next section. 
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3.3 Bundles of multilinear maps over immer- 
sions 



Consider the following natural bundles of k -multilinear mappings: 



Inim 



L'^ (TImm; Tlmm) 



Y 

Imm 



These bundles are isomorphic to the bundles 



L TImm; ] 



Imm 



L TImm; TImm 



Y 

Imm 



where denotes the c°°-completed bornological tensor product of locally con- 
vex vector spaces section 5.7, section 4.29]. Note that iy(TImm; TImm) 
is not isomorphic to T*Imm (g) TImm since the latter bundle corresponds to 
multilinear mappings with finite rank. 

It is worth to write down more explicitly what some of these bundles of 
multilinear mappings are. The tangent space to Imm is given by 

T/Imm = Cf{M,TN) {h G C°^{M,TN) : ttat o /i = /}, 
TImm Cf;?i„,(M, TN) {h e C°°{M, TN) : ttn o h e Imm}. 

Thus T/Imm is the space of vector fields along the immersion /. Now the 
cotangent space to Imm will be described. The symbol (8ic°°(m) nieans that the 
tensor product is taken over the algebra C°°{M). 

Tpmm = L(T/Imm;M) = Cf{M,TN)' = {M)' i§)c'>o^M)Cf{M,T* N) 
T*Imm = L(TImm;M) = C°°(M)' §)c'^^M)CfS,m{M,T* N) 

The bundle L^y^(Tlmm.;M.) is of interest for the definition of a Riemannian 
metric on Imm. (The subscripts sym and ait indicate symmetric and alternat- 
ing multilinear maps, respectively.) Letting (Eis denotes the symmetric tensor 
product and (^s the c°°-completed bornological symmetric tensor product, one 
has 

L^y^,(T/Imm;M) = (T/Imm §s T/Imm)' = {Cf{M,TN) ®s Cf{M,TN))' 
^ {Cf{M,TN (g)s TN))' 
^C°°{M)' ®c^{M)Cf{M,T*N (g)s T*N) 
L,2y^(TImm;M) = C°°(M)' §,c^^M)C^m{M,T* N (g>s T*N) 



3.4. DIFFEOMORPHISM GROUP 
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3.4 DifFeomorphism group 

The following result is taken from section 43.1] with slight simplifications 
due to the compactness of M. 

Theorem. For a smooth compact manifold M the group Diff (A/) of all smooth 
diffeomorphisms of M is an open submanifold ofC°°{M,M). Composition and 
inversion are smooth. The Lie algebra of the smooth infinite dimensional Lie 
group Diff(M) is the convenient vector space X{M) of all smooth vector fields 
on M , equipped with the negative of the usual Lie bracket. DifF(Af ) is a regular 
Lie group in the sense that the right evolution 

evoV : C°°(E,X(Af)) ^ Difr(M) 

as defined in f29l section 38.4] exists and is .smooth. The exponential mapping 

exp : X(Af) ^ Difr(A/) 

is the flow mapping to time 1, and it is .smooth. 

The diffeomorphism group Diff(Af) acts smoothly on C°°{M,N) and its 
subspaces Imm,Immj and Emb by composition from the right. For Imm, the 
action is given by the mapping 

Imm(M, N) x Difr(Af) -> Imm(Af, N), (/, ip) r{f, ip) = r'^(/) = / o ip. 

The tangent prolongation of this group action is given by the mapping 

TImm(M, N) x Diff (M) ^ Tlmm{M, iV), (/i, (/?) ^ Tr^ih) ^ho^p. 

3.5 Riemannian metrics on immersions 

A Riemannian metric G on Imm is a section of the bundle 

L2y,„(TImm;M) 

such that at every / G Imm, G/ is a symmetric positive definite bilinear mapping 

Gf : T/Imm x T/Imm M. 
Each metric is weak in the sense that Gf, seen as a mapping 

Gf : Tyimm — > T^Imm 
is injective. (But it can never be surjective.) 

Assumption. It will always be assumed that the metric G is compatible with 
the action o/ Diff (A/) on Imm(Af, N) in the sense that the group action is given 
by isometrics. 

This means that G = {rf)*G for all ip G Diff(Af), where r'^ denotes the 
right action of (p on Imm that was described in section |3.4[ This condition can 
be spelled out in more details using the definition of r''^ as follows: 



Gf{h, k) = {{r'^YG) (h, k) = Gr^(f) {Tr'f{h),TrV{k)) = Gfo^{h o^,ko^). 
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3.6 Covariant derivative on immersions 



The covariant derivative defined in section |2 . 7l induces a covariant derivative 
over immersions as follows. Let Q be a smooth manifold. Then one identifies 



heC°° (Q, rimm(Af, N)) and X e X{Q) 



with 



h"" eC°"{Qx M, TN) and (X, Om) e X{Q x M). 

As described in section [2.71 one has the covariant derivative 
Vfx,o„)^^eC-(QxM,riV). 

Thus one can define 

Vxh = {^lx.0M)'^''y ^ C°°{Q,Tlmm{M,N)). 



This covariant derivative is torsion-free by section 2.8 formula Q. It respects 
the metric g but in general does not respect G. 

It is helpful to point out some special cases of how this construction can be 
used. The case Q = M will be important to formulate the geodesic equation. The 
expression that will be of interest in the formulation of the geodesic equation 
is Voj/t, which is well-defined when / : M — > Imm is a path of immersions and 
/t : K — >■ TImm is its velocity. 

Another case of interest is Q = Imm. Let h,k,m € X(Imm). Then the 
covariant derivative V„ift. is well-defined and tensorial in to. Requiring V„i to 
respect the grading of the spaces of multilinear maps, to act as a derivation on 
products and to commute with compositions of multilinear maps, one obtains as 
in section [2/7] a covariant derivative Vm acting on all mappings into the natural 
bundles of multilinear mappings over Imm. In particular, VmP and VmG are 
well-defined for 

P e r(L(TImm;TImm)), G G r(L,2y„^(TImm; M)) 

by the usual formulas 

(V,„P)(/i) = V™(P(/i)) - F(V„,/i), 
{WrnG){h, k) = V^{G{h, k)) - G{V^h, k) - G(/i, V„fc). 



3.7 Metric gradients 



The metric gradients H, K € r(L^(TImm; TImm)) are uniquely defined by the 
equation 

(V™G)(/i,fc) = G{K{h,m),k) = G{m,H{h,k)), 



3. 8. GEODESIC EQUATION ON IMMERSIONS 
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where h, k, m are vector fields on Imm and the covariant derivative of the metric 
tensor G is defined as in the previous section. (This is a generahzation of the 
definition used in [32] that allows for a curved ambient space N ^ M".) 

Existence of H, K has to proven case by case for each metric G, usually by 
partial integration. For Sobolev metrics, this will be proven in sections |7.2| and 
[731 

Assumption. Nevertheless it will he assumed for now that the metric gradients 
H, K exist. 



3.8 Geodesic equation on immersions 

Theorem. Given H, K as defined in the previous section and V as defined in 
section \3.6\ the geodesic equation reads as 

y3jt = \Hf{fuft)-Kf{fuft). 

This is the same result as in [351 section 2.4], but in a more general setting. 

Proof. Let / : (— e, e) x [0, 1] x Af — > be a one-parameter family of curves of 
immersions with fixed endpoints. The variational parameter will be denoted by 
s £ (— e,e) and the time-parameter by i S [0, 1]. In the following calculation, 
let Gf denote G composed with /, i.e. 

G/ : M Imm Ll^^iTlumv, M). 

Remember that the covariant derivative on 1mm that has been introduced in 
section [HH] is torsion- free so that one has 



VaJ, - Vajt = Tf.[du d,] + Tor(/i, /,) = 0. 
Thus the first variation of the energy of the curves is 



1 1 



ds^ Gf{ftJt)dt^- iV9Gf){ftJt)+ Gf(yoJtJt)dt 
Jo Jo Jo 

\'^f.G){ftJt)+ r Gf{W9jsJt)dt 



JO 

{VfG){fuft)dt+ [ dt Gf{fsJt)dt 
Jo 

iVf,G){fsJt)dt^ [ Gf{fs,VoJt)dt 



1^ g(/„ ^H{ft, ft) + 0- K(ft, ft) - VoJt)dt. 



If /(O, •, •) is energy-minimizing, then one has at s = that 
1 



^H{ftJt)-KiftJt)-VaJt = 0. □ 
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3.9 Geodesic equation on immersions in terms 
of the momentum 



In the previous section the geodesic equation for the velocity ft has been derived. 
In many apphcations it is more convenient to formulate the geodesic equation 
as an equation for the momentum G{ft, •) G T^Imm. G{ft, •) is an element of 
the smooth cotangent bundle, also called smooth dual, which is given by 

G(rimm) := ]J G'/(TjImm) = ]J {Gf{h, ■) : h e T/Imm} C r*Imm. 

/Gimm /Gimm 

It is strictly smaller than T*Imm since at every / G Imm the metric Gf : 
TfImm — >■ T^Imm is injective but not surjective. It is called smooth since it 
does not contain distributional sections of f*TN, whereas T^Imm does. 

Theorem. The geodesic equation for the momentum p G r*Imm is given by 

P=G{ft,-) 
ya.P^lGf{H{ft,ft),-), 

where H is the metric gradient defined in section \3. 7| and V is the covariant 



derivative action on mappings into r*Imm as defined in section 3.6 



Proof. Let G f denote G composed with the path / : M Imm, i.e. 

G/ : M ^ Imm L^^^(TImm; E). 

Then one has 

Va.P = Va, {Gfift, •)) - (Va,G^)(/i, •) + Gf{VaJt, •) 
= (V/,G)(/i, •) + Gf(^-H{fuft) - K{fu ft), •) 

= Gf{Kifu ft), ■) + Gf [^Hift, ft) - K{ft,ft), •) 



□ 



This equation is equivalent to Hamilton's equation restricted to the smooth 
cotangent bundle: 

f P = G{ft,-) 

|pt = (grad" E){p). 

Here uj denotes the restriction of the canonical symplectic form on T*Imni to 
the smooth cotangent bundle and E is the Hamiltonian 

E : G(TImm) ^ M, E{p) = G-^{p,p) 



which is only defined on the smooth cotangent bundle. 



3.10. CONSERVED MOMENTA 
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3.10 Conserved momenta 



This section describes how a group acting on Imm by isometrics defines a mo- 
mentum mapping that is conserved along geodesies in Imm. It is similar to [71 
section 4]. A more detailed treatment and proofs can be found in |39) . 

Consider an infinite dimensional regular Lie group with Lie algebra q and 
a right action g t-^ oi this group on Imm. Let Imm be endowed with a 
Riemannian metric G. The basic assumption (assumption 3.5 1 is that the action 
is by isometrics: 

G={rSyG, i.e. Gf{h,k)^Gr.if){Tf{r3)h,Tf{rS)k). 

Denote by X(Imm) the set of vector fields on Imm. Then the group action can 
be specified by the fundamental vector field mapping C • ~^ X(Imm), which 
will be a bounded Lie algebra homomorphism. The fundamental vector field 
(^x , X £ Q is the infinitesimal action in the sense: 

C^(/)=at|or°''P(*^)(/). 

The key to the Hamiltonian approach is to write the infinitesimal action as a 
Hamiltonian vector field, i.e. as the w-gradient of some function. This function 
will be called the momentum map. w is a two-form on Tlmm, 

that is obtained as the puUback of the canonical symplectic form on T*Imm via 
the metric 

G : Tlmm — > T*Imm. 
The w-gradient is defined by the relation 

grad" / e X(TImm) , c^(grad'" /, •) = d/, 

where / is a smooth function on Tlmm. Not all functions have an w-gradient 
because 

uj : TTImm — > T*TImm 

is injective, but not surjective. The set of functions that have a smooth lo- 
gradient are denoted by 

C;=°(TImm,M) C C°°(TImm,M). 

The momentum map is defined as 

J : ^ C°°(rimm,M), jx{hf) - Gf{Cx{.f),hf) 

and it is verified that it has the desired properties: Assuming that the metric 
gradients H, K exist (assumption 3.7), it can be proven that 

jx eC^(TImm,M) and grad"(jx) = Cx- 
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Thus the momentum map fits into the following commutative diagram of Lie 
algebras: 



i/0(rimm) ^ C2°('rimm, 




X(TImm, Ld) — "^^-^ H^ (T Imm) 



Here X(TImm, a;) is the space of vector fields on Tlmm whose flow leaves uj fixed. 
All arrows in this diagram are homomorphism of Lie algebras. The sequence at 
the top is exact when it is extended by zeros on the left and right end. 

By Emmy Noether 's theorem, the momentum mapping is constant along any 
geodesic / : M — >■ Imm. Thus for any X G g one has that 

jxift) = Gf{Cx{f), ft) is constant in t. 

Now several group actions on Imm will be considered, and the corresponding 
conserved momenta will be calculated. 



• Consider the smooth right action of the group DifF(M) on Imm(M, N) 
given by composition from the right: 



f^foif for (peDifr(M). 



This action is isometric by assumption, see section 3.5 For X G X{M) 
the fundamental vector field is given by 

a(/) = 9*|o(/oFif) = r/ox 

where Fl^ denotes the flow of X. The reparametrization momentum, for 
any vector field X on M is thus Gf{TfoX, hf). Assuming that the metric 
is reparametrization invariant, it follows that along any geodesic f{t,-), 
the expression Gf{Tf o X, ft) is constant for all X. 

For a flat ambient space N = M" the following group actions can be consider in 
addition: 

• The left action of the Euclidean motion group M" xi SO{n) on Imm(A/, M") 
given by 

f^A + Bf for {A,B) X SO{n). 

The fundamental vector field mapping is 

C(A.x)(/) + for (AX) eM" xso(n). 

The linear momentum is thus G f{A, h), A G M" and if the metric is trans- 
lation invariant, Gf{A, ft) will be constant along geodesies for every A £ 
M". The angular momentum is similarly Gf{X.f,h),X € so{n) and if 
the metric is rotation invariant, then Gf{X.f, ft) will be constant along 
geodesies for each X e so(n). 

• The action of the scaling group of K given by / i— ^ e^ f, with fundamental 
vector field Ca(/) = o,-f- If the metric is scale invariant, then the scaling 
momentum Gf{f,ft) will be constant along geodesies. 



3.11. SHAPE SPACE 
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3.11 Shape space 

Diff(M) acts smoothly on (M, N) and its subsets Imm, Imm/ and Emb by 
composition from the right. Shape space is defined as the orbit space with 
respect to this action. That means that in shape space, two mappings that 
differ only in their parametrization will be regarded the same. 

Theorem. Let M be compact and of dimension < n. Then Imm/(M, iV) is 
the total space of a smooth principal fiber bundle with structure group Diff (M), 
whose base manifold is a Hausdorff smooth Frechet manifold denoted by 

B,j{M, N) = Imm/(M, A^)/Diff(M). 

The same result holds for the open subset Emb(M, A^) C lmnif(M,N). The 
corresponding base space is denoted by 

BeiM,N) = Emb(Af, A^)/Diff(Af). 

However, the space 

B,{M,N) = Imm(A/,iV)/Diff(M) 

is not a smooth manifold, but has singularities of orbifold type: Locally, it looks 
like a finite dimensional orbifold times an infinite dimensional Frechet space. 

The proofs for free and non-free immersions can be found in [T7] and the 
one for embeddings in [29l section 44.1]. 

As with immersions and embeddings, the notation Bij,Bi^Bf, will be used 
when it is clear that M and N are the domain and target of the mappings. 

3.12 Riemannian submersions and geodesies 

The way to induce a Riemannian metric on shape space is to use the concept 
of a Riemannian submersion. This section explains in general terms what a 
Riemannian submersion is and how horizontal geodesies in the top space corre- 
spond nicely to geodesies in the quotient space. The definitions and results of 
this section are taken from [36l section 26]. 

Let TT : — i? be a submersion of smooth manifolds, that is, Ttt : TE TB 
is surjective. Then 

V = V{ti) ker(r7r) C TE 

is called the vertical subbundle. If E carries a Riemannian metric G, then one 
can go on to define the horizontal subbundle as the G-orthogonal complement 
oiV: 

Hor = Hor(7r, G) := ^(Tr)-^ C TE. 

Now any vector X G TE can be decomposed uniquely in vertical and horizontal 
components as 
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This definition extends to the cotangent bundle as follows: An element of T*E 
is called horizontal when it annihilates all vertical vectors, and vertical when it 
annihilates all horizontal vectors. 

In the setting described so far, the mapping 

is an isomorphism of vector spaces for all x G E. If both {E, Ge) and (B, Gb) 
are Riemannian manifolds and if this mapping is an isometry for all x € E, then 
TT will be called a Riemannian submersion. 

Theorem. Consider a Riemannian submersion t: : E ^ B, and let c : [0, 1] — ^ 
E be a geodesic in E. 

1. If c'(t) is horizontal at one t, then it is horizontal at all t. 

2. If c'{t) is horizontal then t: o c is a geodesic in B. 

3. If every curve in B can be lifted to a horizontal curve in E, then there 
is a one-to-one correspondence between curves in B and horizontal curves 
in E. This implies that instead of solving the geodesic equation on B one 
can equivalently solve the equation for horizontal geodesies in E. 

See [311 section 26] for the proof. 



3.13 Riemannian metrics on shape space 

Now the previous chapter is applied to the submersion tt : Imm — )■ Be 

Theorem. Given a DiS{M) -invariant Riemannian metric on Imm, there is a 
unique Riemannian metric on the quotient space Bi such that the quotient map 
TT ; Imm — )• _B, is a Riemannian submersion. 



One also gets a description of the tangent space to shape space: When 
/ G Imm, then T^(^j:-^Bi is isometric to the horizontal bundle at /. The horizontal 
bundle depends on the definition of the metric. For the TJ^-metric, it consists 
of vector fields along / that are everywhere normal to /, see section 5 



Assumption. It will always be assumed that a DiS{M) -invariant Riemannian 
metric on the manifold of immersions is given, and that shape space is endowed 
with the unique Riemannian metric turning the projection into a Riemannian 
submersion. 



3.14 Geodesic equation on shape space 



Theorem |3 . 1 2 1 applied to the Riemannian submersion tt : Imm — !■ Bi yields: 



3.14. GEODESIC EQUATION ON SHAPE SPACE 
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Theorem. Assuming that every curve in Bi can be lifted to a horizontal curve 
in Imm, the geodesic equation on shape space is equivalent to 



ft = /^'- e Hor 



(1) 



lH{fuft)~Kifuft)]''°^ 



where f is a horizontal curve in Imm, where H, K are the metric gradients 
defined in section 3.1 and where V is the covariant derivative defined in sec- 
tion\3.6l 



Proof. Theorem |3. 12| states that the geodesic equation on shape space is equiv- 
alent to the horizontal geodesic equation on Imm which is given by 

r ft = z^'- 

(2) < 1 

[VaJt--i?/(/t,/t)-i^/(/t,/t) 

Clearly ^ implies 0. To prove the converse it remains to show that 

(1 \ vert 

-H{ftJt)-K{fuft)) ■ 

As the following proof shows, this is a consequence of the conservation of the 
momentum along / and of the invariance of the metric under Diff (Af). 

Recall the infinitesimal action of Diff (A/) on Imm(A/, A^). For any X £ 
X(M) it is given by the fundamental vector field 

Cx e X(Imm), Cxif) = 9s|oH/'exp(sX)) = 9,|o(/ o Fl^) = Tf.X. 



Here r is the right action of Diff(A'/) on Imm(Af, A^) defined in section 3.4 
When / : M — )■ Imm is a curve of immersions, one obtains a two-parameter 
family of immersions 

g -.RxR-)- Imm, g{s, t) = r(/(t), exp(sX)) 

that satisfies 

Vojg-ds = VsTg.dt + Tg.[du ds] + Tor(Tg.5t, Tg.d^) 
= Va^r(r'='^P(''-^))/t + + 

since V is torsion-frec. This implies 

Cx{f) is vertical and ft is horizontal by assumption. Thus the momentum 



28 



CHAPTER 3. SHAPE SPACE 



s=0 



mapping Gf[(x{f), ft) is constant and equals zero. Its derivative is 

= dt{Gf{Cx{f)Jt)) 
= iVo,Gf){Cx{f)Jt) +Gf{VeXxif)Jt) +Gf{Cxif),yaJt) 
= (V/,G)(Cx(/),/t) +G/(Va.|„T(r-P(^^))/„/,) +G^(Cx(/),VaJO 
= Gf{Kfif,J,) + V9j„Cx{f)) 

+ G,o.p(.x) / ( Va, r(r'='^p(^^) ) ft , r(r'='^p(^^) ) /t ) [^^ 

^Gf{Kf{ft,ft) + VaJuCx{f)) 

+ ia,|o(^G,..p(..x,^(T(r-P(^^))/„T(r-P(-^))/,) 

- I (V9.G,..p(.x,^)(T(r-P(^^))/i,T(r-P(^^))/4) 
= Gf{Kf{ftJt) + V9jt,Cx{f)) 

+ lds\o{Gfifuft)) ~ ^(Va(/)G)(/,,/t) 

= G/ (^Kfifu ft) + ^aJt + 0- \Hf{fu ft),Cx{f)) 

Any vertical tangent vector to / is of the form Cx{f) for some X G X(M). 
Therefore 

/I \ 
0= (V9jt--i//(/t,/t) + ^/(/t,/t)) • □ 

It will be shown in section |5.9| that curves in Bi can be lifted to horizontal 
curves in Imm for the very general class of Sobolev type metrics. Thus all 
assumptions and conclusions of the theorem hold. 



3.15 Geodesic equation on shape space in terms 
of the momentum 

As in the previous section, theorem |3.12| will be applied to the Riemannian 
submersion tt : Imm — > Bi. But this time, the formulation of the geodesic 



equation in terms of the momentum will be used, sec section 3.9 As will be 
seen in section |5.11[ this is the most convenient formulation of the geodesic 
equation for Sobolev- type metrics. 

Theorem. Assuming that every curve in Bi can be lifted to a horizontal curve 
in Imm, the geodesic equation on shape space is equivalent to the set of equations 

p = G/(/t,-)eHorcT*Imm, 
(Va.p)'^°'-^G/(ff(/„/,),-)'°'-. 

Here f is a curve in Imm, H is the metric gradient defined in section \3.1\ and 
V is the covariant derivative defined in section \3.6[ f is horizontal because p is 
horizontal. 



3.16. INNER VERSUS OUTER METRICS 
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3.16 Inner versus outer metrics 

There are two similar yet different approaches on how to define a Riemannian 
metric on shape space. 

The metrics on shape space presented in this work are induced by metrics on 

Imm(Af, iV). One might call them inner metrics since they arc defined intrin- 
sically to M. Intuitively, these metrics can be seen as describing a deformable 
material that the shape itself is made of. 

In contrast to these metrics, there are also metrics that are induced from 
metrics on DifF(A^) by the same construction of Riemannian submersions. (The 
widely used LDDMM algorithm is based on such a metric.) The difiierential 
operator governing these metrics is defined on all of N, even outside of the shape. 
When the shape is deformed, the surrounding ambient space is deformed with 
it. Intuitively, such metrics can be seen as describing some deformable material 
that the ambient space is made of. Therefore one might call them outer metrics. 

The following diagram illustrates both approaches. Metrics are defined on 
one of the top spaces and induced on the corresponding space below by the 
construction of Riemannian submersions. 



DiS{N) 



Y 



Emb(M, N)^ 



^ Imm(M, N) 



Y 



Be{M,N)' 



^Bi{M,N) 
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Chapter 4 

Variational formulas 



Recall that many operators like 

9 = 1*9, S = Sf, vol(5), V = V^ A = A^ ... 

implicitly depend on the immersion /. In this section their derivative with 
respect to / which is called their first variation will be calculated . These 
formulas will be used to calculate the metric gradients that are needed for the 
geodesic equation. 

This section is based on [T, "S] . Some of the formulas can be found in [T51 
[371 [SU]. The presentation is similar to lOj, and some of the variational formulas 
are the same. 



4.1 Paths of immersions 

All of the differential-geometric concepts introduced in section [2] can be recast 
for a path of immersions instead of a fixed immersion. This allows to study 
variations of immersions. So let / : M — )■ Imm(M, iV) be a path of immersions. 
By convenient calculus [53], / can equivalently be seen as / : M x 7\/ — > iV such 
that f{t, ■) is an immersion for each t. The bundles over M can be replaced by 
bundles over M x M: 

pr* TJM pv^TlM ® f*TN Nor(/) 



MxM MxM MxM 

Here prj denotes the projection prj : M x M — M . The covariant derivative 
y zh is now defined for vector fields Z on M x M and sections h of the above 
bundles. The vector fields (9t, Oa/ ) and (Or, X), where X is a vector field on M, 
are of special importance. In later sections they will be identified with dt and 
X whenever this does not pose any problems. Let 

inst : M ^ M X M, {t, x). 
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Then by property [5] from section |2.7| one has for vector fields X, Y on M 
VxTf{t, ■).Y = VxT(/ o insO o Y = VxTf o Tins* oY 

= VxT/ o (Or, Y) o inSt = VTins, oxT/ o (Or, Y) 

= (V(o„x)r/o(GR,y))oinst. 

This shows that one can recover the static situation at t by using vector fields 
on M X M with vanishing R-component and evaluating at t. 

4.2 Directional derivatives of functions 

The following ways to denote directional derivatives of functions will be used, 
in particular in infinite dimensions. Given a function F{x, y) for instance, 

D[x,h)F will be written as a shorthand for dt\viF{x + t/i, y). 

Here (x, h) in the subscript denotes the tangent vector with foot point x and 
direction h. If F takes values in some linear space, this linear space and its 
tangent space will be identified. 

4.3 Setting for first variations 

In all of this chapter, let / be an immersion and ft £ Tyimm a tangent vector 
to /. The reason for calling the tangent vector /( is that in calculations it will 
often be the derivative of a curve of immersions through /. Using the same 
symbol / for the fixed immersion and for the path of immersions through it, 
one has in fact that 

Du,h)F = dtF{f{t)). 

4.4 Variation of equivariant tensor fields 

Let the smooth mapping F : Imm(M, N) — ?> r(rjM) take values in some space 
of tensor fields over M, or more generally in any natural bundle over M, see 

m- 

Lemma, If F is equivariant with respect to pullbacks by diffeomorphisms of M, 
%. e. 

F{f)^{^*F){f)^^*{F{{^-'rf)) 

for all Lp e Diff(M) and J € Imm(M, A^), then the tangential variation of F is 
its Lie- derivative: 

Du,Tf.f7)F = dt\oF[f o Flf) = dt\^F{{Flfrf) 
= dt\o{Flfy{F{f))^CfT{F{f)). 

This allows us to calculate the tangential variation of the pullback metric 
and the volume density, for example. 



4.5. VARIATION OF THE METRIC 
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4.5 Variation of the metric 

Lemma. The differential of the pullback metric 

Imm ^ T{Slf^T*M), 
_ f -^9^1*9 
is given by 

D^fj^^g = 2 SymgiWfuTf) = ^2g{f,^,S) + 2 Sym V(/7)^ 
^-2g{f,^,S)+Cfrg. 

Here Sym denotes the symmetric part of the tensor field C of type ( 2 ) given 

by 

( Sym(C)) (X, Y) := ^ Y) + C{Y, X)) . 

Proof. Let / : M x Af — > be a path of immersions. Swapping covariant 
derivatives as in section 2.8 formula ^ one gets 



dt{g{X,Y))=dt{g{Tf.X,Tf.Y))^g{yo^Tf.X,TJ.Y)+g{Tf.X,y9jf.Y) 
= gi^xfuTf.Y) + g{Tf.X, Vyft) = (2 ^ymg{V fu Tf)) {X, Y). 
Splitting ft into its normal and tangential part yields 

2 Symg(V/t, T/) = 2 Symg(V/,-^ + VTf.fJ.Tf) 

= -2 Symgift.VTf) + 2 Sym5(V/7, •) 
= -2g(/,^,5) + 2SymV(/7)^ 

Finally the relation 

^(f,Tf.fJ)9 = 2SymV(/7)'' ^ ^fj9 
follows either from the equivariance of g with respect to puUbacks by diffeomor- 



phisms (see section 4.4) or directly from 

{Cxg){Y, Z) = Cx{9{Y, Z)) - g{CxY, Z) - g{Y, CxZ) 

= Vx {9{Y, Z)) - g{VxY - VyX, Z) - g{Y, VxZ - VzX) 

= g{VYX, Z)+g{Y,VzX) = {VyX)\Z) + {V zX)\Y) 

= (VyX^)(Z) + {VzX^){Y) = 2Sym(V(X^))(y,Z). □ 

4.6 Variation of the inverse of the metric 

Lemma. The differential of the inverse of the pullback metric 
( Imm ^ T{L{T*M,TM)), 

I / ^ g-'^iFg)-' 

is given by 

DUM9-' = DuMifg)-' = 2g{ft,g-'Sg-') + Cfr{g-') 
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Proof. 

= 2g-'gif,^, S)g-' - g-\Cjjg)g-' = 2g{J^ , g-^ Sg'^) + C {g-') □ 

4.7 Variation of the volume density 

Lemma. The differential of the volume density 

Imm Vol(M), 



/ ^ vol(5) = vol(/*5) 
is given by 

D(fM^o\ig) = Trf (5(V/t, T/))vol(g) = (divS(/7) - ff(/,^,Tr^(5)))vol(5). 



Proof. Let g{t) € T{SyQT* A4) be any curve of Riemannian metrics. Then 

1 



dtYolig) ^ -Tv{g-\dtg)Yolig). 

This foUows from the formula for vo\{g) in a local oriented chart (u^, . . . u™) on 
M: 



dtvol{g) = dt^Jdet{{g,j)ij) du^ A • • • A du™ 

= , , , Tr(adj (5)9*5) du' A ■ ■ ■ A du"^ 

^ Tr{det{{g^j)^j)g-^dtg) du^ A ■ ■ ■ A du" 



2y^M({g~])~) 
^TTig-\dtg)yol{g) 



Now one can set g — f*g and plug in the formula 

dtg^dt{f*g)^2Symg{Vft,Tf) 
from |4.5| This immediately proves the first formula: 

dtvo\{g) = iTr(.g-i.2Sym5(V/t,T/)) = Tr^ {g{VfuTf)). 

Expanding this further yields the second formula: 

dtY0\{g) = Tr^ {vgifuTf) - g{ft,^Tf)) 

= Ty^ (ygifu Tf) ~ g{ft,S)) = -S/*g{ft,Tf) - 5(/,, Tr3(5)) 

= -V*i{f7f)-g{f,^,TT^{S)) = div(/7) -5(/,^,Tr«(5)). 
Here it has been used that 

VT/ = S and div(/7) = Tr(V/7) = Tr^ ((V/7)^) = -V*((/7)^). 



Note that by |4.4[ the formula for the tangential variation would have followed 
also from the equivariance of the volume form with respect to puUbacks by 
diffeomorphisms. □ 



4.8. VARIATION OF THE COVARIANT DERIVATIVE 
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4.8 Variation of the covariant derivative 

In this section, let V ~ = V-^ ^ be the Levi-Civita covariant derivative 
acting on vector fields on AI. Since any two covariant derivatives on M differ 
by a tensor field, the first variation of V-'' ^ is tensorial. It is given by the tensor 
field D^^fj^)Vf'9 g TmM). 

Lemma. The tensor field D(^f j^-j'V-^ ^ is determined by the following relation 
holding for vector fields X, Y, Z on M: 

9i(Difj,)V){X,Y),Z) = ^{VD(^fj^-jg){X^Y^Z+Y^X®Z-Z®X^Y) 

Proof. The defining formula for the covariant derivative is 

g{VxY, Z) = i [xgiY, Z) + Yg{Z, X) - Zg{X, Y) 

- g{X, [Y, Z\) + [Z, X\) + g[Z, [X, Y\) 
Taking the derivative D{fjt) yields 

(DifMdKyxY, Z) +g{{D^fj^)V)iX, Y), Z) 

- [Y,Z]) + (D^f^gW, [Z,X]) + {D^fj,)g)iZ, [X,Y]) 

Then the result follows by replacing all Lie brackets in the above formula by 
covariant derivatives using [X, Y] — VxY — VyX and by expanding all terms 
of the form X{{D^fj^g){Y, Z)) using 

X{{Du,U)9){Y,Z)) = 

{VxD^jj^)g)iY, Z) + iD^fj^)g){VxY, Z) + {D^fj,)g){Y, VxZ). □ 



4.9 Variation of the Laplacian 

The Laplacian as defined in section |2.11| can be seen as a smooth section of the 
bundle L(TImm; Tlmm) over Imm since for every / G Imm it is a mapping 

A-^ 9 ■ jyimm — >■ Tflmm. 

The right way to define a first variation is to use the covariant derivative defined 
in section [Hm 



Lemma. For A G r(L(TImm; Tlmm)), / G Imm and h G Tyimm one has 
(V„A)(.) = T,(,-..(I.,„„.).,-'V^/.) - V(,,„^, „„^. 

V(R'(f,.Tf)h) - Tr» (^R'(J,,Tf)Vh 
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Proof. Let / be a curve of immersions and h a vector field along /. One has 
A : 1mm L(TImm; Tlmm), Ao/ = A-^*^' :E^Imm^L(TI mm; Tlmm). 
Using property |2.7| 5 one gets 



(VfAm - (Va,(A o f)){h) - Va^Ah ~ AVa^h 
= -Va,TT'{y^h)-AVa^h 

= Tr {g-\D^jj^-,g)g-^S/^h) - Tr^iVa^'h) - A{Va,h). 

The term Tr^(VajV^/i) will be treated further. Let X,Y be vector fields on 
M that are constant in time. When they are seen as vector fields on M x M 
then V^jX = Va^F = 0. Using the formulas from section 2.8 to swap covariant 
derivatives one gets 

(Va,v2/,,)(X,y) = Va,(VxVy/i- Vv^yM 

= VxVa, Vy/i + R^{ft,Tf.X)VYh - Va^v^Yh 
= Vy Va,/i + Vx {R^ifu Tf.Y)h) + RS{ft,Tf.X)VYh 
- '^v^YVd.h- V[a„v^y]/» - R^{fuT f.V xY)h. 

The Lie bracket is 

[dt,VpY]^{D^jj^)V){X,Y) 
since (now without the slight abuse of notation) 

[{duQu), (Oe, vfpY)] = 9,|o TFl^^'f''^ o V^^ o i^^f^OM) 

= (o«,P(/./*)V)(^,n)- 

Therefore 

{Vay^h){XX)^ 

= {V^VaMXX) + V x{R'{fuTf.Y)h) + RS{fuTf.X)VYh 

- ^(DaM^)(x.Y)h - R'ifu Tf.VxY)h 

= {V^WaMX,Y) + {WTf.xR'){ft,Tf.Y)h + R^{Vxft,Tf.Y)h 
+ R^{ft,^xTf.Y)h + Ra{ft,Tf.Y)Wxh + R^{ft,Tf.X)WYh 

- '^[D^,j^,v)ix.Y)h - R^ift, Tf.VxY)h 

= {V^Va,h){X, Y) + {VTf.xR')UuTf.Y)h + R^i^xft, Tf.Y)h 

+ R^{ft, {^Tf){X, Y))h + R^{fu Tf.Y)Wxh + R^{ft,Tf.X)WYh 

~ ^(Du.M^)(X,Y)h 

= {V^VaMX.Y) + Vx{R'{fuTf.Y)h)+R^{fuTf.X)VYh 



4.9. VARIATION OF THE LAPLACIAN 
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Putting together all terms one obtains 

(VfAm = TT{g-\D^fj^)g)g-'V^h) - Tr^ (v {R^{fuTf)h) 

- Tr^ [R^{fuTf)Vh) + VTr.(z5„,,^,v)/i 

= T^{9-\Df^fj^)g)g-'V^h)+V*{R^{JuTf)h) 

- Tr^ (R3{ft,Tf)Vh) + VTr.(D„,,^,v)/i. 



It remains to calculate Tr^(Z)(-yjj) V). Using the variational formula for V from 
section |4.8| one gets for any vector field Z and a g-orthonormal frame Si 



5(Tr«(^(/,/.)V),Z) 

i 

= -{^*iD^f.M9)){Z) - ^Tr^(Vzi?(/j,)5) 
= -(V*(^(/,/.)5))(^) - ^V^Tr^(i?(^,/,)5) 
= -(V*(^(/,A)5) + ^dTr«P(/,/,)5))(^) 
= -5((V*(^(/,/.)5) + ldTr^{D^fj^^g))\ 



Therefore 



Tr^(A/,/*)V) = -{V*iD^fj,)g) + idTr^(Z?(^j,)g))«. □ 
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Chapter 5 

Sobolev-type metrics 



Assumption. Let P be a smooth section of the bundle L(nmm; TImm) over 
Imm such that at every f G Imm the operator 



is an elliptic pseudo differential operator that is positive and symmetric with 
respect to the H'^ -metric on Imm, 



In this section, the geodesic equation on Imm and Bi for the G'^'-metric wih be 
calculated in terms of the operator P and it will be proven that it is well-posed 
under some assumptions. 

This section is based on [8j section 4]. 

5.1 Invariance of P under reparametrizations 

Assumption. It will be assumed that P is invariant under the action of the 
reparametrization group Diff(M) acting on l'nmi{M, N), i.e. 



Pf : T/Imm — >■ Tyimm 




Then P induces a metric on the set of immersions, namely 




M 



g{Pfh, k)vo\{g) for / e Imm, /i, fceT/Imm. 



p = {r^yp 



for all ip e Diff(A/). 



For any / G Imm and (p E Diff(Af ) this means 



Pf - (Tfrn-' o Pfo^ o Tfr 
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Applied to ft, G Tyimm this means 

Pf{h) o (fi = Pfoifih o ip). 

The invariance of P imphes that the induced metric is invariant under 
the action of Diff(Af), too. Therefore it induces a unique metric on Bi as 
explained in section [3. 13| 

5.2 The adjoint of VP 

The following construction is needed to express the metric gradient H which is 
part of the geodesic equation. Hf arises from the metric G/ by differentiating 
it with respect to its foot point / G Imm. Since G is defined via the operator P, 
one also needs to differentiate Pf with respect to its foot point. As for the metric, 
this is accomplished by the covariant derivate. For P e r(L(TImm; TImm)) 
and m G Tlmm one has 

V,„P G r(L(TImm; TImm)) , VP G r{L{T'^lmm; TImm)) . 

See section l376l for more details. 

Assumption. It is assumed that there exists a smooth adjoint 

Adj(VP) G r(i2 (TImm; TImm)) 
of VP in the following sense: 

[ g{{\/,„P)h,k)vo\{g)= [ 5(m,Adj(VP)(ft,fc))vol(5). 

The existence of the adjoint needs to be checked in each specific example, 
usually by partial integration. For the operator P = 1 + AA^, the existence of 
the adjoint will be proven and explicit formulas will be calculated in sections |7.2| 
and El 

Lemma. // the adjoint of VP exists, then its tangential part is determined by 
the invariance of P with respect to reparametrizations: 

Adj(VP)(ft,fc)T = (g(VPft,A;)-5(Vft,Pfc))" 

= grad» g{Ph, k) - {g{Ph, Vfc) + g{Wh, Pk)) " 
for f G Imm, h,k € T^Imm. 

Proof. Let X be a vector field on M. Then 

= Va,|„ {Pf^Fix {h o Flf)) - Pf.pi. (Va,|„ {h o «f )) 
= Va.|o {Pfih) o Flf) - P/(Va,|„(ft o Flf)) 

= Vx{Pf{h))-Pf{Vxh) 
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Therefore one has for m,h,k G Tylmm that 

/ .9(mT,Adj(VP)(/i,fc)^)vol(g)- / 5(r/.mT,Adj(VF)(/i,fc))vol(5) 
= / g{(yTf.m^P)h,k)Y0\{g) 

JM 

- / g(V„T(F/i)-P(V„T/^),fc)vol(g) 

JM 

{g{V^TPh,k)-g{V„,Th,Pk))vo\{g) 

M 



M 



g(mT, (g(VP/i, fc) - giVh, Pfc))") vol(g). 



5.3 Metric gradients 



As explained in section |3.8[ the geodesic equation can be expressed in terms of 
the metric gradients H and K. These gradients wiU be computed now. 

Lemma. // Adj(VP) exists, then also H and K exist and are given by 



Kf{h,m) = F-i((V„,P)/i + Tr9 {g{\/m,Tf)).Phj 
Hf{h,k) = p-i(Adj(VP)(/i,fc)^ -P/.(5(P/i,Vfc) +g(V/i,PA:))" 
-g{Ph,k).Tr<^{S)] 



Proof. For vector fields m, h, k on Imm one has 

(V„,G^)(/i, k) = / g{Ph, fc)vol(5) 



(1) 



M 

5(P(V™/i),fc)vol(5)- / g{Ph,\/^k)vo\{g) 

M JM 

= I D^f^,^)g{Ph,k)vo\{g)+ I g(P/i, fc)D(y,„,)Vol(.g) 

JM JM 

g(P(V™/i),fc)vol(5)- / g{Ph,\/^k)vo\{g) 

M JM 

= [ g{V,niPh),k)vo\{g)+ [ g{Ph,V „,k)Yo\{g) 

JM JM 



g{Ph, k)D(f^„-,)Vol{g) 

M 

g{PiVmh),k)vo\{g)- [ ^(P/i, V™fc)vol(g) 

M JM 

g{{\7rnP)h,k)vo\{g)+ [ g{Ph,k)D(f^^)Yo\{g) 

M JM 

One immediately gets the if- gradient by plugging in the variational formula |4.7| 
for the volume form: 

Kf{h, m) = PJ^ ({VmP)h + Tr9 (g(Vm, T/)) .Ph 
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To calculate the i?-gradient, one rewrites equation ([T]) using the definition of 
the adjoint: 

{V,nG''){h,k)^ f 5(m,Adj(VP)(/i,fc))vol(.g)+ / g{Ph,k)D(j,^)Vo\{g). 

Now the second summand is treated further using again the variational formula 
of the volume density from section |4.7[ 



{Ph, fc)i?(/,,„)V0l(5) = / g{Ph, k) Tr^ (^(Vm, T/))vol(5) 

M J hi 

= / g{Ph,k)TvB {Vg{m,Tf) -g{m,VTJ))vo\{g) 

J M 

g{Ph,k)(^- V*g(m,T/) -g(m,Tr9(S'))^vol(g) 



= -/ g"i{Vg{Ph,k),g{m,Tf))vo\{g)- g [P h, k)g{m, Ti^ iS))vo\{g) 
= / g{m, -Tf. grad^ g{Ph, k) - g{Ph, k) Tr3(S'))vol(.g) 

J M 

Collecting terms one gets that 
GP{Hf{h,k),m)^{VmG^){h,k) 

= / g(m,Adj(VP)(/i,fc)-r/.grad^5(P;i,fc)-5(P/i,fc)Tr3(5))vol(<?) 
Jm 

Thus the i?-gradient is given by 

Hf{h, k) = P-I ( Adj(VP)(/i, k) - Tf. grad^ g{Ph, k) - g{Ph, k).TT^S)^ 

The highest order term grad^ g{Ph, k) cancels out when taking into account the 
formula for the tangential part of the adjoint from section [5?2} 



Hfih, k) = P-I ( Adj(VP)(/i, k)^ - Tf.{giPh, Vfc) + g(V/i, Pfc))" 

-g{Ph,k).TTa{S)). □ 



5.4 Geodesic equation on immersions 

The geodesic equation for a general metric on Imm(M, N) has been calculated 
in section [TSl and reads as 

Plugging in the formulas for H, K derived in the last section yields the following 
theorem. 

Theorem. The geodesic equation for a Sobolev-type metric on immersions 
is given by 

^ojt =^P-i ( Adj(VP)(/„ /O-" - 2.Tf.g{Pfu V0 - g{Pft, ft).TT^S) 
- P-'((.Vf,P)ft + Tr^ {giVfuTf))Pft 
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5.5 Geodesic equation on immersions in terms 
of the momentum 

The geo desic equation in terms of the momentum has been calculated in sec- 
tion 3.9 for a general metric on immersions. For a Sobolev-type metric , the 



momentum G^{ft, ■) takes the form 

Pft(g) vol(g) : M r*Imm 

since all other parts of the metric (namely the integral and g) are constant and 
can be neglected. 

Theorem. The geodesic equation written in terms of the momentum for a 
Sobolev-type metric G^ on Imm is given by: 

p = Pft® vol(.g) 

Va,p = ^ ( Adj(VP)(/„ /O^ - 2r/.5(P/„ V/,)« - /t) Tr«(^)) ® vol(g) 



5.6 Well-posedness of the geodesic equation 

It will be proven that the geodesic equation for a Sobolev-type metric G^ on 
Imm is well-posed under some assumptions on P. It will also be shown that 
(tt, exp) is a diffeomorphism from a neighbourhood of the zero section in TImm 
to a neighbourhood of the diagonal in Imm x Imm. 

First, Sobolev sections of vector bundles are introduced. More information 
can be found in [44^ and in ^2] ■ Let E ^ M he a. vector bundle over a compact 
Riemannian manifold (M, g). E is given a fiber Riemannian metric g^ and a 
compatible covariant derivative on r{E) is chosen. Let V*^ be the Levi- 
Civita covariant derivative on T{TM) for a fixed background metric g on M . 
Then the Sobolev space H^{E) is the Hilbert space completion of the space of 
smooth sections V{E) in the Sobolev norm 




The Sobolev space does not depend on the choices of g,g ,V ; the resulting 
norms are equivalent, see f331. The following results hold (see pT]): 

Sobolev lemma. If k > dim(M)/2 then the identy on T{E) extends to a 
injective bounded linear mapping H^'^^^E) — ^ Cp{E) where C^^E) carries the 
supremum norm of all derivatives up to order p. 

Module property of Sobolev spaces. If k > dim(A/)/2 then pointwise eval- 
uation H^{L{E^E)) X II^[E) — > II^{E) is bounded bilinear. Likewise all other 
pointwise contraction operations are multilinear bounded operations. 



44 



CHAPTER 5. SOBOLEV-TYPE METRICS 



The following notation shall be used: 

rimm(A/, N) = Ci^^{M, TN) := {s G C°°(M, TN) : ttn o s e Inim(M, N)}, 
TflTJim{M,N) = Cf{M,TN) := {s e C^{M,TN) : ttat o s = /} ^ T{f*TN). 

The smooth Sobolev manifolds (for k > dim(Af)/2 + 1) 

Imni''+^ (Af, N) C Imm'' (M, N) , Q Imm'' (M, N) = Inim(Af, N) 

k 

shall also be considered. They are constructed from the Sobolev completions in 
each canonical chart separately and then glued together, always with respect to 
the the background metrics. See [21] for a detailed treatment in a more general 
situation (Af does not need to be compact there). 

It is assumed that the operator P satisfies the following properties, for k > 
diin(M) _|_ ^ (^QjjiQ of the assumptions have already been stated earlier.) 

Assumption 1. P e r(L(nmm; TImm)) is smooth and invariant under the 
action o/Diff(A/). (See section \5.1\ for the definition of invariance.) 

Assumption 2. For each / e Imm(A/, N) the operator 

Pf : T{f*TN) T{f*TN) 

is an elliptic pseudo- differential operator of order 2p for p > of classical type 
which is positive and symmetric with respect to the H'^ -metric on Imm, 

H^f{h, k) = g{h, fc)vol(g) for f G Imm and h,k £ Tylmm. 

Since Pf is elliptic, it is unbounded selfadjoint on the Hilbert completion of 
Tflmm with respect to see [HI theorem 26.2]. Furthermore Pf extends to 
a bounded injective (since it is positive) linear operator 

Pf : H^+^P(f*TN) ^ H^{f*TN) 

which is also surjective (since it is Fredholm as an elliptic operator, with van- 
ishing index as a selfadjoint operator). 

Assumption 3. P extends to a smooth section of the smooth Sobolev bundle 
L «mm^+^p iM, TN) , H'^^^^, {M, TN) \ ) 

I 

lYa.m*'+'^P{M,N) 

with fiber 

L{H''+^P{f*TN),H''{f*TN)) over f e lmm''+^P{M, N) 

such that 

Pf : H''+^'P{f*TN) H''{f*TN) 
is invertible for each f G Imm''^^^(Af, A^). 
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By the implicit function theorem on Banach spaces, f ^-^ Pj.- ^ is then a 
smooth section of the smooth Sobolev bundle 



L{H^^^, (M, A^) |i„„,.+..(M,Ar) , (M, N)) 



Imm'''+2P(M,iV) 



Moreover, ^ is also an elliptic pseudo-differential operator of order — 2p. 
Assumption 4. T/ie normal part of the adjoint f i— )■ Adj(VP)^ from sec 



tion (5.2) extends to a smooth section of the smooth Sobolev bundle 
L%UH'+''if*TN);HHf*TN)) 



Y 

Imm'^+^P^^^^) 

Assumption 5. All mappings Pju, PJ^u, and Adj(VP)^(M, are linear 
pseudo differential operators in u and v of order 2p, ~2p, and 2p, respectively. 
As mappings in the footpoint f , viewed locally in a trivialisation of the bun- 
dle Cf^^{M,TN) 1tom\{M,N), they are nonlinear, and it is assumed that 
they are a composition of operators of the following type: (a) Local operators 
of order <2p, i.e., A{f){x) ^ A{x,V^p f{x),i/^P-^ f{x), . . . ,V fix), f{x)). (b) 
Linear pseudo-differential operators. 

These properties hold for the operator 1 + A/S.^ considered in section [7) 

Theorem. Let p > 1 and k > dim(Af)/2 + 1, and let P satisfy assumptions 
1-5. 



Then the initial value problem for the geodesic equation (5.4) has unique 
local solutions in the Sobolev manifold Imm'^'*'^'' of H^^'^p -immersions. The 
solutions depend smoothly on t and on the initial conditions /(O, . ) and ft{0, . ). 
The domain of existence (in t) is uniform in k and thus this also holds in 
Imm(M, N). 

Moreover, in each Sobolev completion Inim''"'"^^, the Riemannian exponen- 
tial mapping exp^ exists and it smooth on a -open neighborhood of the 
zero section in the tangent bundle, and (tt, exp^) is a diffeomorphism from a 
(smaller) H^"-open neigbourhood of the zero section to an H'^"-open neighbor- 
hood of the diagonal in Imm'^"'^^^ x Imm*^^^'', where fcg is the smallest integer 
> dim(A/)/2 + 2p. All these neighborhoods are uniform in k > dim(M)/2, and 
thus both properties of the exponential mapping continue to hold in Imm(M, iV). 



This proof is partly an adaptation of [35", section 4.3]. It works in three steps: 
First, the geodesic equation is formulated as the flow equation of a smooth 
vector field on a Sobolev completion of Tlmm. Thus one gets local existence 
and uniqueness of solutions. Second, it is shown that the time-interval where a 
solution exists does not depend on the order of the Sobolev space of immersions. 
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Thus one gets solutions on the intersection of all Sobolev spaces, which is the 
space of smooth immersions. Third, a general argument proves the claims about 
the exponential map. 



Proof. The geodesic equation is considered as the flow equation of a smooth 
{C°°) autonomous vector field X on 

Locally near a fixed smooth immersion /q e Imm(M, N) this looks like JJ'^^'^p x 
H^U^TN) where C/'^+^p is i/'^-open in Imm'=+2P(M, TV) C H''+^p{M,N). It 
suffices to prove the theorem in each open subset f/'^+^P x H^{fQTN). 



On this subset one can write X = {Xi,X2) as follows (using (5.4)) 
ft^Pi'u^:X,{f,u) 



ut = ^Adi{VP)iP-'u,p-'u)^ 



(5) -Tf.g{u,VPy'uy - ^g{u,p-'u).Trf'^Sf) 



For (f,u) e U'^+^P X H'' one has Adi{W P) f{Py^u, Py'^u) e . When Y e 
H'^ then also Y^ e since g = f *g G i/'^+^P-i c H^^'^. Similarly, when 
a € H'' then also a"^ e H''. & 77^+2^-2 ^ ^/c_ -pj^^g ^^^.^^^ ^^^.^^ 

investigation using (1) ~ (3) shows that the right hand side of (5) is smooth in 
(/, u) € U^+'^P X H'^ with values in H^+'^p x . Thus by the theory of smooth 
ode's on Banach spaces, the flow Fl'^' exists on [/'=+2p x and is smooth in 
t and the initial conditions for fixed k > '^™^(-'^'^) ^ 

Consider C°° initial conditions /o = /(O, ) and mq = Pj^ftiO, ) = 
u(0, ) for the flow equation (5) in U°° x r(/o TiV). Suppose the trajectory 
Flj (/o, Wo) of X through these initial conditions in U^+'^p x maximally exists 
for t G (-Ofe, 5fc), and the trajectory FlJ^+^/o, uo) in [/'=+i+2p x H''+^ maximally 
exists for t G {—ak+i,bk+i) with ak+i < and 6^+1 < b^, say. By uniqueness 
of solutions one has Fl("'"^(/o, uq) = Fl^(/o, uq) for t e (— a/£+i,fofc+i). Now V is 
applied to both equations (5): 

(V/)t=V/t=VXi(/,u) 
iVu)t^Vut = VX-2if,u) 

It is claimed that the highest derivatives of / and u appear only linearly in 
VX,{f,u) for i = 1,2, i.e. 

VX,(/, u) - «)(V2P+i/) + X,,2(/, u)(V2f+iu) + X,,3(/, u) 

where all Xij{f,u){v) and Xi^3{f,u) ~ 1,2) are smooth in all variables, of 
highest order 2p in / and u, linear and algebraic (i.e., of order 0) in v. This 
claim follows from assumption [5] (a) For a local operator we can apply the 
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chain rule: The derivative of order 2p + 1 of / appears only linearly, (b) For 
a linear pseudo differential operator A of order k the commutator [V, A] is a 
pseudo-differential operator of order k again. 

Then one writes V^^+^Z — V^^/ and V^^+^m — V^^'u for the highest deriva- 
tives only. The last system now becomes 

/t = Xi4(/,U)(V2P/) + Xi,2(/,^/)(V2Pu) +Xi,3(/,«) 
Ut = X2,l(/,u)(V2f/) +X2,2(/,U)(V2?'U) +X2^3(/,U) 

which is inhomogeneous bounded linear in {f,u) G {7*=+2p x H'' with coeffi- 
cients bounded linear operators on H^+'^p and , respectively. These coeffi- 
cients are C°° functions of (/, u) £ jj^+'^p x which are already known on 
the interval {—ak,bk)- This equation therefore has a solution {f{t, ),u(t, )) 
for all t for which the coefficients exists, thus for all t G {ak,bk)- The limit 
\imty^bk+i{f{ti )) exists in U^+'^p x and by continuity it equals 

(V/, Vu) in f/'^+^P X for some i > fefe+i. Thus the i7''+^-flow was not max- 
imal and can be continued. So (— a/j+i, fe^+i) = (— afe,6fe). Iterating this proce- 
dure one concludes that the flow of X exists in nm>fc C/™^^^ x i/™ = Imm x C°° . 

It remains to check the properties of the Riemannian exponential mapping 
exp^. It is given by exp^(/i) — c(l) where c{t) is the geodesic emanating 
from value / with initial velocity h. The properties claimed follow from local 
existence and uniqueness of solutions to the geodesic equation on each space 
Imm'^^^^(M, A^) and from the form of the geodesic equation fu = ^f{ft,ft) 
when it is written down in a chart using the Christoffel symbols, namely linearity 
in /(( and bilinearity in ft- See for example [3BJ 22.6 and 22.7,] for a detailed 
proof in terms of the spray vector field S{f, h) = (/, h; h,r f{h, h)) which works 
on each TImm''+^P(M, N) = H^^'^^^^^{M, TN) without any change in notation. 

So one checks this on the largest of these spaces Imm*"'" (M, A^) (i.e. with the 
smallest k). Since the spray on Imm'^'^ restricts to the spray on each Imm*^^^^, 
the exponential mapping exp^ and the inverse (7r,exp^)~^ on Imm*^" restrict 
to the corresponding mappings on each Imm*'"''^^. Thus the neighborhoods of 
existence are uniform in k. □ 



5.7 Momentum mappings 

Recall that by assumption, the operator P is invariant under the action of the 
reparametrization group Diff (M) . Therefore the induced metric is invariant 
under this group action, too. According to [Mi section 2.5] one gets: 

Theorem. The reparametrization momentum, which is the momentum map- 
ping corresponding to the action o/DifF(Af) on Imm(il/, iV), is conserved along 
any geodesic f in Imm(Af, iV); 



\fXeX{M) : [ g{Tf.X,Pf,)vo\{g) 
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or equivalently 

g{{PftV)Yol{g) e r(r*M®M vol(M)) 

is constant along f. 



5.8 Horizontal bundle 



The splitting of Tlmm into horizontal and v ertica l subspaces will be calcu- 
lated for Sobolev-type metrics . See section 3.12 for the general theory. By 
definition, a tangent vector h to f £ Imm(A/, N) is horizontal if and only if 
it is G^-perpendicular to the Diff (Af )-orbits. This is the case if and only if 
g{Pfh{x),T^f.Xx) = at every point x e M. Therefore the horizontal bundle 
at the point / equals 

{h e T/Imm : Pfh{x)±T^f{T^M) for all x € M} ^ {h : [PfhY = O}. 

Note that the horizontal bundle consists of vector fields that are normal to / 
when P = Id, i.e. for the 7?°-metric on Imm. 

Let us work out the G^-decomposition of h into vertical and horizontal parts. 
This decomposition is written as 

h = T/./i™'' + h^°'. 

Then 

Pj:h = Pf{Tf.h:'''') + Pfh^°' with {PfhY = {Pf{Tf.h^''')y + 0. 

Thus one considers the operators 

Pj : X(A/) ^ X{M), Pj(X) = {Pf{Tf.X)f, 

Pf^i_ : X(M) ^ r(Nor(/)) C G°°(Af,TiV), = {Pf{Tf.X))^. 

The operator Pj is unbounded, positive and symmetric on the Hilbcrt comple- 
tion of T/Imm with respect to the iJ'^-metric since one has 



/ g{PjX,Y)yoUg)= [ g{Tf.PjX,Tf.Y)yoUg) 

JM Jm 

g{Pf{Tf.X),Tf.Y)vo\{g) 



M 



g{PjY,X)vo\{g), 



f g{PjX,X)Yol{g)^ f g{Pf{Tf.X),Tf.X)vol{g) >0 ii X ^ 0. 

JM JM 

Let a^^ and a^f denote the principal symbols of Pf and Pj , respectively. Take 
any x G AI and £, G T*M \ {0}. Then {^) is symmetric, positive definite on 
(Ty(2,-)iV, g). This means that one has for any h,k G Tf(^^-^N that 

g{aPf iOh, k) = g{h, o^' (Ok) , gia""' {OK h) > for h ^ 0. 
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The principal symbols a^f and a^f are related by 

g{<jPhOX,Y) ^giT/.a^'f {OX,Tf.Y) ^g{a'''{OTf.X,Tf.Y), 

where X, F e T^M. Thus a^f (^) is symmetric, positive definite on (T^M^g). 
Therefore Pj is again elliptic, thus it is selfadjoint, so its index (as operator 
jjk+2p _^ jjk-^ vanishes. It is injective (since positive), hence it is bijective and 
thus invertible. Thus it has been proven: 

Lemma. The decomposition of h ^ Tlmm into its vertical and horizontal com- 
ponents is given by 

5.9 Horizontal curves 

To establish the one-to-one correspondence between curves in shape space and 
horizontal curves in Imm that has been described in theorem |3.12[ one needs 
the following property: 

Lemma. For any smooth path f in Imm(Af , N) there exists a smooth path ip 
in DifF(M) with <^(0, . ) — IcLm depending smoothly on f such that the path f 
given by f{t,x) = f(t,ip{t,x)) is horizontal: 

g{Pjidtf),Tf.TM)^0. 

Thus any path in shape space can be lifted to a horizontal path of immersions. 

The basic idea is to write the path ip as the integral curve of a time dependent 
vector field. This method is called the Moser-trick (see [171 Section 2.5]). 

Proof. Since P is invariant, one has (r^)* P ~ P or Pfo^{u o ip) ~ {Pfu) o ip 
for ip e DifF(M). In the following fop will denote the map f{t,ip{t,x)), etc. 
One looks for tp as the integral curve of a time dependent vector field ^(t, x) on 
M, given hy (pt = o ip. The following expression must vanish for all x € M 
and e T^M: 

Q^g(Pfo^{dt{fo^)){x),T{fop).X,) 
= g{Pfo^{{dtf) o p + Tf.idt^)) (x), T{f o 
= g{{{PfidJ)) + PfiTf.O) {^ix)),Tf o Tp,.X,)j 
Since Tip is surjective, Tip.X exhausts the tangent space T^p(^x^M, and one has 

{{Pf{dtf))+Pf{Tf.O){p^{x)) ± f. 
This holds for all x G Af , and by the surjectivity of ip, one also has that 
{{Pf{dJ)) + Pf{Tf.O){x) ± f 
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at all X £ M. This means that the tangential part [Pf{dtf) + Pf{Tf.£)) 
vanishes. Using the time dependent vector field 

i = -{pJ)-\{Pfd,fY) 

and its flow ip achieves this. □ 



5.10 Geodesic equation on shape space 



By the previous section and theorem |3.12[ geodesies in Bi correspond exactly to 
horizontal geodesies in Imm. The equations for horizontal geodesies in the space 



of immersions have been written down in section 3.14 Here they are specialized 
to Sobolev-type metrics: 

Theorem. The geodesic equation on shape space for a Sobolev-type metric 
is equivalent to the set of equations 



ft 



f'r e Hor, 



kd\{^P){fuft)^ - g{Pfu ft).'TT^{S) 



P-H{{ynP)ft)^ - Tr« {g{Wft,Tf))Pft 



where f is a horizontal path of immersions. 



These equations are not handable very well since taking the horizontal part 
of a vector to Imm involves inverting an elliptic pseudo-differential operator, see 
section 5.8 However, the formulation in the next section is much better. 



5.11 Geodesic equation on shape space in terms 
of the momentum 

The geodesic equation in terms of the momentum has been derived in sec- 



tion 3.15 for a general metric on shape space. Now it is specialized to Sobolev- 



type metrics using the formula for the i?-gradient from section 5.3 



As in section 5.5 the momentum G^{ft, •) is identified with Pft (g) vol(g). 
By definition, the momentum is horizontal if it annihilates all vertical vectors. 
This is the case if and only if Pft is normal to /. Thus the splitting of the 
momentum in horizontal and vertical parts is given by 

Pft ® volCg) = (Pft)^ (g> vol(.g) + Tf.{Pft)^ (g, Yo\{g). 

This is much simpler than the splitting of the velocity in horizontal and vertical 
parts where a pseudo-differential operator has to be inverted, see section |5.8[ 
Thus the following version of the geodesic equation on shape space is the easiest 
to solve. 
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Theorem. The geodesic equation on shape space is equivalent to the set of 
equations for a path of immersions f : 

p = Pft®yo\{g), Pft = iPft)^, 

(Va.p)'^" = \ (Adi{VPM Jt)^ ~g{PftJt).Tr^S)) ®yo\{g). 

The equation for geodesies on Imm without the horizontahty condition is 
\/o,p - ^( Adj(VP)(/i, ft)^ - 2Tf.g{Pfu ^ft)^ - g(Pfu ft)- ^^^{S)) ® vol(g), 



see section 5.5 It has been proven in section 3.14 that the vertical part of this 
equation is satisfied automaticaUy when the geodesic is horizontal. Nevertheless 
this will be checked by hand because the proof is much simpler here than in the 
general case. 

If ft is horizontal then by definition Pft is normal to /. Thus one has for 
any X e X(M) that 

g{{Va,Pft)^.X)=g{Ve,PfuTf.X) = Q-g{PfuVeJf.X) 
= -giPfuVxft) = -g{g{Pft,Vft)^X). 

Thus 

(Ve.p)™''' = ((Va,P/0 ® vol(5) + Pft ® i5(/,/,)vol(g))™'-* 

= Tf.{\Io,PftV ® vol(5) + Tf.{Pft)^ (g> D^fj^)V0\{g) 
= -Tf.g{Pft,VftY®yo\{g)+Q, 

which is exactly the vertical part of the geodesic equation. 
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Chapter 6 

Geodesic distance on shape 
space 

It came as a big surprise when it was discovered in |38| that the Sobolev metric 
of order zero induces vanishing geodesic distance on shape space Bi. It wiU 
be shown that this problem can be overcome by using higher order Sobolev 
metrics. The proof of this result is based on bounding the G'^-length of a path 
from below by its area swept out. The main result is in section [6. 6[ 

This section is based on [51 section 5]. The same ideas can also be found in 
pm section 2.4], [7, section 7] and [37, section 3]. 



6.1 Geodesic distance on shape space 

Geodesic distance on Bi is given by 

dist^-p(Fo,Fi) = infi^i.(F), 

F 

where the infimum is taken over all F : [0, 1] — > Bi with F{0) — Fq and F{1) — 
Fi . L^^p is the length of paths in Bi given by 

ig*^ [F) = ^ ^GP{Ft,Ft)dt for F : [0, 1] ^ B,. 
Letting tt : Imm — )■ Bi denote the projection, one has 

L^'pin o /) = L^--(/) ^ ^GP{fuft)dt 

when / : [0, 1] — >■ Imm is horizontal. In the following sections, conditions on the 
metric G^ ensuring that dist^'p separates points in Bi will be developed. 
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6.2 Vanishing geodesic distance 

Theorem. The distance dist^o induced by the Soholev metric of order zero 
vanishes. Indeed it is possible to connect any two distinct shapes by a path of 
arbitrarily short length. 

This resuh was first estabhshed by Michor and Mumford for the case of 
planar curves in [35]. Here a more general version from [33 is quoted. 



Proof. Take a path f{t,x) in Imm(M, N ) fro m /o to /i and make it horizontal 



by the same method that was used in 5.9 Horizontality for the iJ"-metric 



simply means g{ft,Tf) = 0. This forces a reparametrization on /i. 

Let a : M [0, 1] be a surjective Morse function whose singular values are 
all contained in the set : < fc < 2A^} for some integer N. We shall use 
integers n below and we shall use only multiples of N. 

Then the level sets Mr := {x E M : a{x) — r} are of Lebesgue measure 
0. We shall also need the slices M^^^ra ■— {x S M : ri < a{x) < r2}. Since 
M is compact there exists a constant C such that the following estimate holds 
uniformly in t: 

f vol(/(i, Yg) < C{r2 - ri) / vo\{f{t, )*!)) 

J Mr, „ JM 



Let f{t,x) — f{ip{t,a{x)),x) where f : [0,1] x [0,1] — > [0,1] is given as in 
I 3.10 by 



(p{t,a) = 



'2t{2na~2k), Q < i < 1/2, f < a < ^ 

2t{2k + 2- 2na), < i < 1/2, ^ <a< ^ 

2i - 1 + 2(1 - t){2na - 2k), 1/2 < < < 1, |^ < a < ^ 

^2i - 1 + 2(1 - t){2k + 2- 2na), 1/2 <t< 1, ^ <a< 



See figure |6.1| for an illustration of the construction. 

Then we get T f = ip^.da.ff + Tf and ft — y^t-ft where 

-4nt ( Ana - 4fc 

-Ant I 4fc + 4 - Ana 

" ^ +An{l -t) ' '^^ ^ p- Ana + Ak 
-An{l-t) [-(2-4na + 4fc) 

We use horizontality g{ft,Tf) = to determine f^ = ft + Tf{X) where X G 
TM satisfies = g{ft + T/(X), T/(0) for aU ^ e TM. We also use 



da{0 = 5(grad^ a, = ^(^/(grad^ a),Tf{0) 
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Figure 6.1: Plot of the function (p. Each zig-zagged hne corresponds to fit,-) 
for some fixed values of t, namely t — j^, ^g, . . . , j^. The figure is taken from 

m- 



and get 

= '^t.(/?a-ff(grad^a,OII/«lP+ 

+ ^^<7(grad^ a, X).5(grad« a, OWftf + g{Tf{X), T/(0) 
= (ipt.^c. + fl-gigre^d^ a, f ^(grad^ a, + g{X, 

This implies that X = X grad^ a for a function A and in fact we get 



- ^ f, Tf fgrad^ a) 



and 



ftf 



l + ^l\\da\\j,JM^ 
From Tf — ipa.da.ft + Tf and g{ft,Tf) = we get for the volume form 
vol(fg) = .^l + ^l\\da\\l\\ft\\^ vol(5). 
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For the horizontal length we get 

Jo ^Jm ' 



1 



M^l + <^2||rf^||2||;^||2 



vol(g) ) ^ dt 



V^^-^^2fc 2fe+i AA + (4nt)'IM«ll^ll/dl' 

2n • 2ra * 

+ / ^i±i^i22lMLvoi(.))y.i+ 

•^^2fc+l 2fe+2 Jl + (4nt)2||da||2||/,||2 V 
2n ' 2n V 

^i£^i2^±aMLvoi(.H 

. V^oVM^fe 2fc+i Jl + (4n(l-t))2||rfa||2||/,||2 

2n' 2n V 

•^^2fc+l 2fc+2 ./l + (4n(l-f))2||da||2||/,||2 ! } 

Let £ > 0. The function (t, .x) a(a;)), is uniformly bounded. On 

2k+i the function Ana — 4k has values in [0, 2]. Choose disjoint geodesic 

2ti ' 2ri 

balls centered at the finitely many singular values of the Morse function a of total 
fif-volume < s. Restricted to the union Mging of these balls the integral above is 
0(l)e. So we have to estimate the integrals on the complement M = Af\Msing 
where the function ||(iQ;||g is uniformly bounded from below by ry > 0. 

Let us estimate one of the sums above. We use the fact that the singular 
points of the Morse function a lie all on the boundaries of the sets M2k 2fc+i 

2n' 2n 

so that we can transform the integrals as follows: 

n-l 



E/ ^i^2£^aMLvoi(,) 

T^Jm^^^ yi + (4nt)2||da||2||/,||2 



2n • 2n 
n-l „2fc+l 



y f (4nr-4fc)^||/t|p vol(Crff) 

We estimate this sum of integrals: Consider first the set of all [t, r, x) G Mr 
such that \ft{ip{t,r),x)\ < e. There we estimate by 

0(l).n.l6n2.£2.(^3/3)|r=l/2n ^ 

On the complementary set where \ ft{<f{t,r), x)\ > s we estimate by 
0(l).n.l6n2.— V(^V3)i;-y'" = 0{^^) 
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which goes to if n is large enough. The other sums of integrals can be estimated 
similarly, thus L)^™{f) goes to for n — >■ oo. It is clear that one can approximate 
</3 by a smooth function without changing the estimates essentially. □ 



6.3 Area swept out 

For a path of immersions / seen as a mapping / : [0, 1] x M — > iV one has 
(area swept out by /) = / vol(/(-, ■)*!)) = / / ||/t^|| vol{g)dt. 

J[0,l]xM Jo JM 

6.4 Area swept out bound 

Lemma. Let he a Sobolev type metric that is at least as strong as the i/"- 
metric, i.e. there is a constant Ci > such that 



WHgp > Ci \Mh« J ^ TImm. 

Then one has the area swept out bound for any path of immersions f : 
Ci {area swept out by f) < max y^Vol(/(t)) .LJ^?'(/). 



The proof is an adaptation of the one given in [3 section 7.3] for almost 
local metrics. 



Proof. 



L'^'\.f) = t ll/tllcp dt > Ci ll/J^o dt 
Jo Jo 

>Ci f'Wft^Hodt^C^ r ( I ||/,^f vol(5))^ 
Jo Jo ^Jm ' 

>C, f \ f vol(5))"' / l.||/,^||vol(.g)di 
-'0 ^Jm ' Jm 



dt 



>Cimin( / vol(5)) ' / vol(/(., O^s) 

* ^Jm ' J[0,l]xM 

— CA max / vol(g) j ^ (area swept out by /) □ 
^ * Jm ' ' 
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6.5 Lipschitz continuity of VVol 

Lemma. Let be a Sobolev type metric that is at least as strong as the H^- 
metric, i.e. there is a constant C2 > such that 



\\h\\a^ > C2 = J g{{l + A)h,h)vo\{g) for all h e Tlmm. 



Then the mapping 
is Lipschitz continuous, i.e. for all Fq and Fi in Bi one has: 



/Vol : (Bi,dist^i.) ^ M>o 



^Yo\{Fi) - VVol(i^o) < ;r^dist^'p(Fo,Fi). 



For the case of planar curves, this has been proven in [321 section 4.7]. 
Proof. 

dt\o\ = j^^ (div^(/7) - g{ft,TT^{S))yo\{g) 

= 0+ / 5(/t,V*T/)vol(5)= / (g?®5)(V/t,r/)vol(5) 

JM JM 



Thus 



< y j^^ l|V/t||^o^^vol(5)^ j^^ ||T/||^o^^vol(5) 

< Vvol< -L||/J^, Vvol. 

(-'2 



By integration one gets 



^ 1 „„„ 1 



Now the infimum over all paths / : [0, 1] — > Imm with 7r(/(0)) = Fq and 
7r(/(l)) = Fi is taken. □ 



6.6 Non-vanishing geodesic distance 

Using the estimates proven above and the fact that the area swept out separates 
points at least on Be, one gets the following result: 
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Theorem. The Sobolev type metric induces non-vanishing geodesic distance 
on Be if it is stronger or as strong as the -metric, i.e. if there is a constant 
C > such that 



> C \\h\\jj^ " / ^((^ + for all h e Tlmm. 
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Chapter 7 



Sobolev metrics induced by 
the Laplace operator 



The results on non-vanishing geodesic distance from the previous section lead 
us to consider operators P that are induced by the Laplacian operator: 

P 1 + AAP, P e r(L(rimm; TImm)) 



for a constant A > 0. (See section 2.11| for the definition of the Laplacian 



that is used in this work.) At every / G Inim, Pf is a positive, selfadjoint 
and bijective operator of order 2p acting on T^Imm = r{f*TN). Note that 
A depends smoothly on the immersion / via the puUback-metric f*g, so that 
the same is true of P. P is invariant under the action of the reparametrization 
group Diff(Af). It induces the Sobolev metric 

G^{h,k)^ ( g{Pf{h),k)yo\{g)^ ( 5((l + A(A/*T)^fc)vol(rg). 

JM JM ^ ' 

When ^ = 1 we write i/^ := G^^^" . 

In this section we will calculate explicitly for P = 1 + ^A^ the geodesic 
equation and conserved momenta that have been deduced in section [5] for a 
general operator P. The hardest part will be the partial integration needed for 
the adjoint of VP. As a result we will get explicit formulas that are ready to 
be implemented numerically. 

This section is based on [5J section 6]. 



7.1 Other choices for P 

Other choices for P are the operator P = 1 + yl(V*)''V^ corresponding to the 
metric 

G;(/i,fc)= / (5(/i,fc) + A5(VP/i,Vffc))vol(g), 

JM 
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and other operators that differ only in lower order terms. Since these operators 
all have the same principal symbol, they induce equivalent metrics on each 
tangent space T/Imm. It would be interesting to know if the induced geodesic 
distances on Bi are equivalent as well. 



7.2 Adjoint of VP 

To find a formula for the geodesic equation one has to calculate the adjoint 
of VP, see section |5.4| The following calculations at the same time show the 
existence of the adjoint. It has been shown in section [5.2| that the invariance 
of the operator P with respect to reparametrizations determines the tangential 
part of the adjoint: 

Adj(VP)(/i,fc))^ = grad^ g(P/i, fc) - {g{Ph,Vk) +g{Vh, Pk)f . 

It remains to calculate its normal part using the variational formulas from sec- 
tion H 

In the following calculations there will be terms of the form TT{g~^ sig~^ S2) , 
where si,S2 are two-forms on M. When the two-forms are seen as mappings 
TM T*M, they can be composed with g^^ : T* M -> TM. Thus the expres- 
sion under the trace is a mapping TAI — >■ TM to which the trace can be applied. 
When one of the two-forms is vector valued, the same tensor components as be- 
fore are contracted. For example when h G r{f*TN) then S2 = V^/i is a 
two-form on M with values in f*TN. Then in the expression Ti{g~^ .si.g~^ .82) 
only TM and T*M components are contracted, whereas the f*TN component 
remains unaffected. 



5(m^,Adj(VP)(/i,fc))vol(5)= / g((V„,^P)/i, fc)vol(5) 

M JM 



= Aj2f 5((V^.A)Af-^'i/z,A'fc)vol(ff) 

= ^E / g(Tr(g-\D^f^^^)g.g-'V^AP-^~'h),A^k)Yo\ig) 
4=0 -^A^ ^ ^ ^ ^ 

-Ay^ qhf , s,AP-'-i/i,A*fcVol(ff) 

^A/^ (v*(I?(,,„^,9) + i<iTr«(D,,,„,,ff)) ; 

+ AV / g(v*P>i^,P/)Af-*-i/i,A'fc)vol(5) 
4=0 -^A^ ^ ^ 

-AV / g(Tr«(P^(m^,r/)VAP-*-i/i),A^fc)vol(5) 
/ Tr(g-i.i^(/,,„.)5.5-ig(V2AP-'-i/i,A^fc))vol(5) 



7.2. ADJOINT OF VP 
p-i 



M 

1 



-AJ2 (5?®3)(vA^'"*-i/i,(V*i?(;,„,.)g)®A*fc)vol(5) 

P-1 j. , 

-AY, / (5?§5g)(vAf-*-i/i,-rfTr3(Z?(/,„.)5)® A'fc)vol(g) 
p-i ^ 

+ / (5?®g)ta™^,r/)AJ'"^-i/i,VA^A:)vol(5) 

p-i „ 

-AV / g(Tr3(i?^(m^,T/)VAf-^-i/i),A^fc)vol(5) 

Using the following symmetry property of the curvature tensor (see |36l 24.4. 
g{B?{X, Y)Z, U) = -9{R^{Y, X)Z, U) - -g{R!^{Z, U)Y, X) 

yields: 

g(TO^,Adj(VP)(/i,fc))vol(.g) = 
p-i . 

/ .92(^(/.m-)ff,5(V2AP-^-i;.,A^fc))vol(5) 
p-1 „ 

-A^ / g?(5(VAP-^-i/i,A^fc),V*i^(/,,„.)5)vol(g) 

-^E / 5?(5(VAP-^-i/z,A^fc) ivTr^(i?(/,„.)5))vol(g) 
i=o •'^ ^ 
p-i „ 

+ ^E / ff(Tr^ (^^(A^"'"^/i,VA*fc)T/),m-L)vol(5) 
i=o -^^^ 
p-i f 

-AV / gfxr^ (i?^(VAP-'-i/i,A*fc)T/),m-L)vol(5) 
p-i . 

^E / 52(^(/.m-)ff,5(V2AP''-i/i,A^fc))vol(5) 
p-1 „ 

-^E / ff2(v5(VAf— i;j,A'fc),7^(/,„.)5)vol(g) 

- ^ E / (v*5(VA''-^-i/i, A^fc)) TrS(i?(^,^.)5)vol(g) 
i=o "^^-^ 
p-1 „ 

+ ^E / ^(Tr^ (^''(A^"'"^/i,VA*fc)T/),m-L)vol(5) 

P-1 . 

-AV / gfxr^ (i?^(VAP-*-i/i, A*fc)T/),m-L)vol(5) 
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= AJ2 / 52°(^(/.r„-)3,5(V'Af-'-i/i,A'A;))vol(ff) 

p-i „ 

-^E / 52(5(V'AP-'-i/i,A'fc),D(/,^.)g)vol(5) 
-^E / 52"(5(VAf-^-i/i,VA^fc),Z)(/,„x)5)vol(ff) 

- 4 E / (v*5(VAf-'-i/i, A^fc)) Tr^(Z)(/,„.)5)vol(ff) 

+ AV / fff Tr^ (i?9(Af-'-i/i,VA'fc)T/),m-L)vol(5) 

-^V / g(TrO{R^{VAP-'-'h,A'k)Tf),m^)Yol{g) 

= -^5^ / 52°(^(/,m-)5,5(VAf-'-i/i,VA'A;))vol(ff) 
. p— 1 „ 

- 2 E / (v*5(VAf-'-i/i, A'fc)) Tr«(I)(/,„.)5)vol(5) 
i=o -'^ 

p-i „ 

+ AV / ^(Trs (i?9(Af-'-i/i,VA'fc)T/),m-L)vol(5) 

P-l r. 

-^V / 5(Tr»(i?^(VAf-^-i/i,A'/c)T/),m^)vol(5) 

p-i . 

= -^E / 52(-2.5(m^,5),5(VAP-'-i/i,VA'fc))vol(ff) 

- I E / (v*5(VAf-'-i/i, A^fc)) TVS ( _ 2.g(rn^, S))yo\{g) 

p-i „ 

+ ^E / ^(Trs (i?9(AP-*-i/i,VA*fc)r/),m-L)vol(5) 

p-i „ 

-^E / ^(Trs (i?^(VAJ'-*-i/i,A*fc)r/),m-L)vol(5) 

„ p-i 
= / g(m^,2Ay^Tr{g-'Sg-^g{VAP-'-^h,VA'k))) 

„ p-i 
+ / 5(m^,^E(^*^(VAf-^-ift,A*fc))Tr'?(5))vol(5) 
•^A^ ^ i=o ^ 

p-i „ 

+ ^E / 5(Tr9(i?'^(A^'-^"'/i,VA^fc)T/),m^)vol(5) 
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^ V / g(TT^RS{VAP-'-'h,A^k)Tf),', 



m Uo\{g). 



From this, one can read off the normal part of the adjoint. Thus one gets: 

Lemma. The adjoint ofVP defined in section \57^ for the operator P = l+^A^ 
is 

p-i 

Adj(VP)(/i,/c) = 2A^Tr(5-i5'5-i5(VAP-*-i/i,VA^fc)) 



1=0 



p-i 



i=0 

p-i 

+ A^TrS {RB{AP^'-^h,\/A'k)Tf) 

i=Q 
p-1 

- A^Tr^ {R^{VAP-'-^h,A'k)Tf) 

4 = 

+ Tf. [ grad^ g[Ph, k) - {g{Ph, Vk) + g{Vh, Pk)) ^ 



7.3 Geodesic equations and conserved momen- 
tum 



The shortest and most convenient formulation of the geodesic equation is in 



terms of the momentum p = (1 + AAP)ft (g) vol(g), see sections 5.5 and 5.11 

Theorem. The geodesic equation on Imm(Af, A^) for the -metric with P 
1 + AA^ is given by: 

f p={l + AAP)ft<E>vo\{g), 
/ p-i 

Ve,P - U E {g-'Sg-'g{W{AP—'ft), VA7t)) 
V j=0 



i=0 

2AETr« (i?nA^-^~Vt,VA7t)r/) 



1=0 



\g{PIu ft) Tr^(5) - Tf.g{PfuVftA <E> vo\{g). 



This equation is well-posed by theorem 5.6 since all conditions are satisfied. For 



the special case of plane curves, this agrees with the geodesic equation calculated 
in f3^ section 4.2]. 

P = 1 + ^A'' and consequently G^ are invariant under the action of the 



reparametrization group Diff(Af). According to sections 3.10 and 5.7 one gets: 
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Theorem. The momentum mapping for the action ofT)iS{M) on Imm(A{f, N) 



+ AAP)ft) ) ^ vol(5) e r(r*M vo1(M)) 
is constant along any geodesic f in Inim(M, A^). 



The horizontal geodesic equation for a general metric on Imm has been 
derived in section [3.15[ In section [5.11| it has been shown that this equation takes 
a very simple form. Now it is possible to write down this equation specifically 
for the operator P = 1 + AA^: 

Theorem. The geodesic equation on shape space for the Sobolev-metric 
with P ^ I + AA^ is equivalent to the set of equations 

p = Pft®Y0\(g), Pft = {Pft)^, 
(Va.p)'^- = A^Tr(g-i5g-i5(VAf— i/t,VA7t)) 



A 



p-i 



i=0 

p-i 



Tr^ {RS{AP-'-'fu \/A'ft)Tf) 



i=0 



-g{PftJt).Ti'S{S)]®yo\{g), 



where f is a path of immersions. For the special case of plane curves, this agrees 
with the geodesic equation calculated in ]39l section 4-S]- 



Chapter 8 

Surfaces in n-space 



This section is about the special case where the ambient space N is M". The 
flatness of M" leads to a simplification of the geodesic equation, and the Eu- 
clidean motion group acting on M" induces additional conserved quantities. The 
vector space structure of K" allows to define a Frechet metric. This metric will 
be compared to Sobolev metrics. Finally in section [83] the space of concentric 
hyper-spheres in K" is briefly investigated. 

Most of the material presented here can also be found in [5]. 



8.1 Geodesic equation 

The covariant derivative on is but the usual derivative. Therefore the 
covariant derivatives V^j/t and V^^p in the geodesic equation can be replaced 
by ftt and pt, respectively. (Note that Imm(M, E") is an open subset of the 
Frechet vector space C°°(M, K").) Also, the curvature terms disappear because 
M" is flat. Any of the formulations of the geodesic equation presented so far can 
thus be adapted to the case N = R^. 

We want to show how the geodesic equation simplifies further under the 
additional assumptions that dim(M) = dim(7V) — 1 and that M is orientable. 
Then it is possible define a unit vector field v to M. The condition that ft is 
horizontal then simplifies to Pft = a.v for a S C°°{M). The geodesic equation 
can then be written as an equation for a. However, the equation is slightly 
simpler when it is written as an equation for a.Yo\{g). In practise, vol(g) can 
be treated as a function on M because one can identify vol(g) with its density 
with respect to du^ A ... A du"^^, where (m^, . . . , is a chart on M. Thus 

multiplication by vol{g) does not pose a problem. 

Theorem. The geodesic equation for a Sobolev-type metric on shape space 
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i?i(M, M") with dim(Af) = ?i — 1 is equivalent to the set of equations 
Pft - 

dt{a.yo\ig)) = ig( Adj(VP)(/t, /t), J^) - \giPfu fMTi^iS),,^), 
where f is a path in Inim(7\f, M") and a is a time- dependent function on M. 



Proof. Applying g{-Tv) to the geodesic equation 5.11 on shape space in terms 
of the momentum one gets 

dt{a.Yo\{g)) = dt g{Pft ® vol(g), v) 

= gi^d^ {Pft ® vol(5)) , + g{Pft ® Y0\{g)yQ^v) 

= \g{Ad]{\/P){fuft),v) ~ \g{Pfuft)g{'T^'{S),v)+Q. □ 

Let us spell this equation out in even more details for the iJ^-metric. This 
is the case of interest for the numerical examples in section |10[ 

Theorem. The geodesic equation on shape space Bi{M,W^) for the Sobolev- 
metric with P = 1 + AA is equivalent to the set of equations 

Pft = a.iy 

9,(a.vol(g)) = {Ag",{s,giVfuVft)) ^ ( + ^ II V/t||^c^^ )) vol(g), 

where f is a path of immersions, a is a time- dependent function on M , s = 
g{S,iy) € r(r2A/) is the shape operator, L = g~^s G T{TlM) is the Weingarten 
mapping, and Tr(L) is the mean curvature. 



Proof. The fastest way to get to this equation is to apply g{-, v) to the geodesic 
equation on Imm from section |7.3[ This yields 



9t(a.vol(5)) - (AT^{g-\s.g~'g{Vft.yft)) + | (V*5(V/t, /t)) Tr(L) 



[g-\s.g~^g[yft,yh)} + ^ 
i5(i^/t,/t).Tr(L))vol(5) 



= (Agl{s,-gy f,y f,)) - ^Tr« (^(V^/,, /,)) Tr(L) 

- ^ Tr'^ (5(V/t, V/0) Tr(i) - ^^((l + AA)/,, /,) Tr(L)) vol(g) 

= [Agl{s,gy fty ft)) - ^ Tr^ {g{Vft,Vft)) Tr(L) 
-^5(/t,/t)-Tr(L))vol(.9). 

Notice that the second order derivatives of ft have canceled out. □ 
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8.2 Additional conserved momenta 

If P is invariant under the action of the EucHdean motion group M" x SO(n), 
then also the metric is in invariant under this group action and one gets 



additional conserved quantities as described in section 3.10 



Theorem. For an operator P that is invariant under the action of the Euclidean 
motion group R" xi SO(n), the linear momentum 



Pftvo\{g) e {Wr 

M 



and the angular momentum 



VX€so(n) : / g{X.f,Pft)vo\{g) 

J M 



or equivalently I (/A Pft)vo\{g) e /\^M" = So{n)* 
J M 

are constant along any geodesic f in Imm(Af, M"). The operator P = 1 + AA^ 
satisfies this property. 

These momenta have also been calculated in [§1 section 4.3]. 



8.3 Frechet distance and Finsler metric 

This section can also be found in fS^, section 5.6]. 

The Frechet distance on shape space i3i(M, M") is defined as 
dist^'(^^o,i^i)= inf ll/o-AILo., 

Jo Jl 

where the infimum is taken over all /o,/i with 7r(/o) = Fo,Tr{fi) = Fi. As 
before, tt denotes the projection tt : Imm — )• Bi. Fixing /o and /i, one has 



distS(^(/o),^(/i)) =inf||/oo(^-/i|| 



where the infimum is taken over all (p E Diff (M). The Frechet distance is related 
to the Finsler metric 



G°° : rimm(Af,M") ^ M, \\h-^\ 



Lemma. The pathlength distance induced by the Finsler metric G°° provides 
an upper bound for the Frechet distance: 



dist^' (Fo , Fi ) < dist^'^ (Fo , Fi ) = inf / 1 1 /, | 

J Jo 



GOO dt, 



where the infimum is taken over all paths 

/: [0,l]^Imm(Af,M") with ^(/(O)) = Fq, 7r(/(l)) = Fi. 
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Proof. Since any path / can be reparametrized such that ft is normal to /, one 
has 



J Jo J Jo 



where the infimum is taken over the same class of paths / as described above. 
Therefore 



distS(Fo,Fi) = inf 11/(1) - /(0)|1^„ = inf 



ftdt 



= inf 
/ 



Jo 



dt = dist^^{Fo,Fi) 



<inf / ll/tlLo 



dt 



□ 



It is claimed in [3^, theorem 13] that c?oo — distcoo. However, the proof 
given there only works on the vector space C°°(M, K") and not on i3j(M, K"). 
The reason is that convex combinations of immersions are used in the proof, 
but that the space of immersions is not convex. 



8.4 Sobolev versus Frechet distance 



This section can also be found in [8, section 5.7]. 

It is a desirable property of any distance on shape space to be stronger 
than the Frechet distance. Otherwise, singular points of a shape could move 
arbitrarily far away without increasing the distance much. 

As the following result shows, Sobolev metrics of low order do not have this 
property. The author and the authors of ^ believe that they have this property 
when the order is high enough, but were not able to prove this. 

Lemma. Let be a metric on _Bi(M, M") that is weaker than or at least as 
weak as a Sobolev -metric with p < '^'™(-'^^) -|- j 



\\h\\ap <C\\h\\H, = J g{{l + AP)h,h)vol{g) for all h € Tlmm. 
Then the Frechet distance can not be bounded by the -distance. 



Proof. It is sufficient to prove the claim for _P = 1 + A^. Let /o be a fixed 
immersion of M into M", and let /i be a translation of fo by a vector h of 
length £. It will be shown that the 7?^-distance between 7r(/o) and 7r(/i) is 
bounded by a constant 2L that does not depend on £, where tt denotes the 
projection of Imm onto Bi. Then it follows that the iJP-distance can not be 
bounded from below by the Frechet distance, and this proves the claim. 

For small tq, one calculates the iJP-length of the following path of immer- 
sions: First scale fo by a factor tq, then translate it by h, and then scale it 
again until it has reached /i. The following calculation shows that under the 
assumption p < m/2 + 1 the immersion /o can be scaled down to zero in finite 
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iyP-pathlength L. Let r : [0, 1] — > [0, 1] be a function of time with r(0) — 1 and 
r(l) = 0. 

L?Z^{r-fo) = ^^^g(r,.(l + (A('--/o)*9)p)(/o),ri./o)vol((r./o)*5)rfi 

= ^^^r2.g((l + -^(A/o*9)f)(/o),/o)r"vol(/*g)dt 

= ^^^^((l + ^(A/o-9)f)(/o),/o)r"vol(/*g)dr =: L 

The last integral converges if < ^1, which holds by assumption. Scaling 

down to ro > needs even less effort. So one sees that the length of the shrinking 
and growing part of the path is bounded by 2L. 

The length of the translation is simply ^■\/f™Vol(/o) = 0(r™/^) since the 
Laplacian of the constant vector field vanishes. Therefore 

distgf (7r(/o), 7r(/i)) < distf„l(/o, /i) < 2L. □ 



8.5 Concentric spheres 

This section can also be found in [8, section 6.6]. 

For a Sobolev type metric that is invariant under the action of the SO{n) 
on M", the set of hyper-spheres in M" with common center is a totally geodesic 
subspace of 5^(5'""^, M"). The reason is that it is the fixed point set of the group 
SO{n) acting on Bi isometrically. (One also needs uniqueness of solutions to 
the geodesic equation to prove that the concentric spheres are totally geodesic.) 
This section mainly deals with the case P = 1 + A^". 

First we want to determine under what conditions the set of concentric 
spheres is geodesically complete under the G^-metric. 

Lemma. The space of concentric spheres is complete with respect to the 
metric with P = 1 + AAp iff p > (n + l)/2. 



Proof. The space is complete if and only if it is impossible to scale a sphere 
down to zero or up to infinity in finite G^ path-length. So let / be a path 
of concentric spheres. It is uniquely described by its radius r. Its velocity is 
/( = rt.v, where v designates the unit normal vector field. One has 

g{g-\S,u)=:L = -lldTM, Tr(L'=) - (-l)'^^^!, Vol - r"-i p^J^. 



Keep in mind that r and rt are constant functions on the sphere, so that all 
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derivatives of them vanish. Therefore 

Aiy ^ V*(Vi/) = V*{-Tf.L) ^ Tr^ (y{Tf.L) 
= Tr^ (v(r/).L) + Tr^ (r/.(VL)) 
= Tr(L2).^ + Tr^ (T/.V(-i Mtm)) - "^-i^ + 

and 

{n-l)P 



Pft = il + AAP){rt.:y)^rt. [1 + A 



From this it is clear that the path / is horizontal Therefore its length as a path 
in Bi is the same as its length as a path in Imm. One calculates its length as 



in the proof of 8.4 



, , , An - l)P\ n.7r"/2 



2P J r(n/2 + 1) 



The integral diverges for ri — !• oo since the integrand is greater than r'"~^)/^. 
It diverges for ro — ^ iff (n — 1 — 2p)/2 < —1, which is equivalent to p > 
(n+l)/2. □ 

The geodesic equation within the space of concentric spheres reduces to an 
ODE for the radius that can be read off the geodesic equation in section |7.3[ 

'n~l p.A.{n-lY 



2/' 



2r r(r^P + A{ji - 1)P) 



Chapter 9 

Diffeomorphism groups 



For M — N the space Emb(M, M) equals the diffeomorphism group of M. An 
operator P g r(L(TEmb; TEmb)) that is invariant under reparametrizations 
induces a right-invariant Riemannian metric on this space. Thus one gets the 
geodesic equation for right-invariant Sobolev metrics on diffeomorphism groups 
and well-posedness of this equation. To the authors knowledge, well-posedness 
has so far only been shown for the special case M = 5'"'^ in [18] and for the 



special case of Sobolev order one metrics in [23]. Theorem 5.6 establishes this 



result for arbitrary compact M and Sobolev metrics of arbitrary order. 

In the decomposition of a vector h e T/Emb into its tangential and normal 
components h = Tf.hJ -\- h-^, the normal part h-^ vanishes. By the naturality 
of the covariant derivative, S = VT/ vanishes as well. Thus on gets: 

Theorem. The geodesic equation on Diff(M) is given by 

f p = Pft<E)voUg), 

\ Va,P = -Tf.g{Pfu V ftf ® vol(ff). 



Note that this equation is not right-trivialized, in contrast to the equation 
given in [5l|371|35], for example. 
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Chapter 10 

Numerical results 



This section can also be found in [5^ section 7]. 

It is of great interest for shape comparison to solve the boundary value prob- 
lem for geodesies in shape space. When the boundary value problem can be 
solved, then any shape can be encoded as the initial momentum of a geodesic 
starting at a fixed reference shape. Since the initial momenta all lie in the same 
vector space, this also opens the way to statistics on shape space. 

There are two approaches to solving the boundary value problem. In [71 the 
first approach of minimizing horizontal path energy over the set of curves in Imm 
connecting two fixed boundary shapes has been taken. This has been done for 
several almost local metrics. For these metrics it is straightforward to calculate 
the horizontal energy because the horizontal bundle equals the normal bundle. 
However, in the case of Sobolev type metrics the horizontal energy involves the 



inverse of a differential operator (see section 5.81, which makes this approach 
much harder. 

The second approach is the method of geodesic shooting. This method is 
based on iteratively solving the initial value problem while suitably adapting 
the initial conditions. The theoretical requirements of existence of solutions to 
the geodesic equation and smooth dependence on initial conditions are met for 



Sobolev type metrics, see section 5.6 This makes geodesic shooting a promising 
approach. 

The first step towards this aim is to numerically solve the initial value prob- 
lem for geodesies, at least for the H^-metric and the case of surfaces in M.^, and 
that is what will be presented in this work. 

The geodesic equation on shape space is equivalent to the horizontal geodesic 
equation on the space of immersions. For the case of surfaces in M.^, it takes 
the form given in section |8.1[ This equation can be conveniently set up using 
the DifferentialGeometry package incorporated in the computer algebra system 



Maple as demonstrated in figure 10.1 (The equations that have actually been 
solved were simplified by multiplying intermediate terms with suitable powers 
of -\/voI(g) , but for the sake of clearness this has not been included in the Maple 
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code in figure 10.1 ) 



Unfortunately, Maple (as of version 14) is not able to solve PDEs with more 
than one space variable numerically. Thus the equations were translated into 
Mathematica. The PDE was solved using the method of lines. Spatial discretiza- 
tion was done using an equidistant grid, and spatial derivatives were replaced 
by finite differences. The time-derivative ft appears implicitly in the equation 
Pfift) = a-i^, and this remains so when the operator Pf is replaced by finite 
differences. 

The solver that has been used is the Implicit Differential- Algebraic (IDA) 
solver that is part of the SUNDIALS suite and is integrated into Mathematica. 
IDA uses backward differentiation of order 1 to 5 with variable step widths. 
Order 5 is standard and has also been used here. At each time step, the new 
value of ft is computed using some previous values of /, and then the new value 
of / is calculated from the equation Pfift) = The dependence on / in 

this equation is of course highly nonlinear. A Newton method is used to solve 
it. This operation is quite costly and has to be done at every step, which is 
a main disadvantage of backward differentiation algorithms. Explicit methods 
are probably much better adapted to the problem. The implementation of 
an explicit solver is ongoing work of Martin Bauer, Martins Bruveris and the 
author. The results are too preliminary to be included in this thesis. 

In the examples that follow, / at time zero is a square [0, tt] x [0, tt] flatly 
embedded in M.^. This is a manifold with boundary, but it can be seen as a 
part of a bigger closed manifold. Zero boundary conditions are used for both / 
and a. It remains to specify an initial condition for a. As a first example, let 
us assume that a at time zero equals sin(x) sin(?/), where x, y are the Euclidean 
coordinates on the square. The resulting geodesic is depicted in figure |10.2[ 
In the absence of a closed-form solution of the geodesic equation, one way to 
check if the solution is correct is to see if the energy G^{ft,ft) is conserved. 
Figure |10.3| shows this for the geodesic from figure |10.2| with various space and 
time discretizations. 



A more complicated example of a geodesic is shown in figure 10.4 and 10.5 
There, the initial velocity was chosen to be a smoothened version of a black 
and white picture of the letter A. The initial momentum was computed from it 
using a discrete Fourier transform. 



Finally, it is shown in figure 10.6 that self-intersections of the surface can 
occur. This is not due to a numerical error but part of the theory, and can be 
an advantage or a disadvantage depending on the application. 
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1 with(Dif f erentialGeometry) ;with(Tensor) ;with(Tools) ; 

DGsetup ( [u , v] , [x , y , z] , E) ; 
3 f := evalDG(f l(t,u,v)*D_x+f2(t,u,v)*D_y+f3(t,u,v)*D_z) ; 

G := evalDG(dx &t dx + dy &t dy + dz &t dz) ; 
5 Gaiimia_vrt := femult Connection (dx &t D_x &t du) ; 

Tf := CovariantDerivative(f ,Gamina_vrt) ; 
7 g:=ContractIndices(G &t Tf &t Tf , [ [1 , 3] , [2 , 5] ] ) ; 

g_inv:=InverseMetric(g) ; 
9 Gaiimia_bas : =Christof f el(g) ; 

det : =Hook ( [D_u , D_u] , g) *Hook ( [D_v , D_v] , g) -Hook ( [D_u , D_v] , g) ~2 ; 
11 ft := evalDG(diff (f l(t,u,v) ,t)*D_x+diff (f2(t,u,v) ,t)*D_y 

+diff (f3(t,u,v) ,t)*D_z) ; 
13 ft := convert (ft .DGtensor) ; 

S := CovariantDerivative(Tf ,Gamma_vrt ,Gairana_bas) ; 
15 cross :=evalDG((dy &w dz) &t D_x + (dz few dx) &t D_y + (dx few dy) &t D_z) ; 

N:=Hook( [ContractIndices(Tf &t D_u, [ [2,3] ] ) , 
17 ContractIndices(Tf &t D_v, [ [2,3] ])], cross) ; 

nu:=convert(N/sqrt(Hook([N,N] ,G)) ,DGvector) ; 
19 s := ContractIndices(G &t S &t nu, [ [1 ,3] , [2 , 6] ] ) ; 

L := Contractlndices (g_inv &t s,[[2,3]]); 
21 Gftft := Contractlndices (G &t ft &t f t , [ [1 ,3] , [2 ,4] ] ) ; 

Cft := CovariantDerivative (f t ,Gaimna_vrt) ; 
23 CCft := CovariantDerivative (Cf t ,Gamma_vrt ,Gaiimia_bas) ; 

Dft := Contractlndices (-g_inv &t CCf t , [ [1 ,4] , [2 , 5] ] ) ; 
25 Pft := evalDG(ft+A*Dft) ; 

GCftCft := Contractlndices (Cft &t Cft &t G, [ [1 ,5] , [3, 6] ] ) ; 
27 gGCftCft := Contractlndices (g_inv &t GCftCft, [[2,3]] ) ; 

TrLgGCftCft := Contractlndices (L &t gGCf tCf t , [ [1 ,4] , [2 , 3] ] ) ; 
29 TrgGCftCft := Contractlndices (gGCf tCft ,[ [1 , 2] ]) ; 

TrL := Contractlndices (L, [ [1 , 2] ]) ; 
31 # b(t,u,v) := a(t,u,v)*sqrt(det) ; 

eql := Contractlndices (Pft &t dx, [ [1 , 2] ] ) *sqrt (det) = 
33 b(t,u,v)*ContractIndices(nu &t dx,[[l,2]]); 

eq2 := Contractlndices (Pft &t dy , [ [1 , 2] ] ) *sqrt (det) = 
35 b(t,u,v)*ContractIndices(nu &t dy,[[l,2]]); 

eq3 := Contractlndices (Pft &t dz , [ [1 , 2] ] ) *sqrt (det) = 
37 b(t,u,v)*ContractIndices(nu &t dz,[[l,2]]); 

eq4 := diff (b(t,u,v) ,t) = 
39 (A*TrLgGCftCft - TrL/2* (Gf tf t+A*TrgGCf tCf t) ) *sqrt (det) ; 

Figure 10.1: Maple source code to set up the geodesic equation. 



Figure 10.2: Geodesic where a bump is formed out a flat plane. The initial 
momentum is a = sin(a;) sin(?/). Time increases linearly from left to right. The 
final time is ^ = 5. 
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Fig ure 1 0.3: Conservation of the energy G^{ft,ft) along the geodesic in fig- 
The true value of G^{ft,ft) is 



10.2 



0.822467 for ^ = 1. The 



4(1+2A) 

maximum time step used in blue and green is 0.1. For purple and cyan it is 
0.05. The number of grid points used in blue and cyan is 100 times 100. For 
green and purple it is 200 times 200. 



Figure 10.4: Letter A forming along a geodesic path. Time increases linearly 
from left to right. The final time is t — 0.8. Top and bottom row are different 
views of the same geodesic. 
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Figure 10.5: Initial velocity /t(0, •) and momentum a(0, •) of the geodesic in 
figure [TO^ Both are shown first from above, then from the side. 



^'^ ^^^^^ 



Figure 10.6; A self-intersection forming along a geodesic. Time increases linearly 
from left to right. 
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Zusammenfassung 



ZuSAMMENFASSUNG. Geometrische Figuren spiclcn cine wichtige Rolle in vie- 
len Bereichen der Wissenschaft, Technik und Medizin. Mathematisch kann eine 
Figur als unparametrisierte, immersive Teil-Mannigfaltigkeit modelliert werden, 
was auch dor dicscr Arbeit zugrunde liegende FigurenbegrifF ist. Wenn der Raum 
der Figuren mit einer Riemannschen Metrik ausgestattet wird, eroffnet sich 
die Welt der Riemannschen Differential-Geometrie mit Geodaten, Gradienten- 
Fliissen und Kriimmung. Leider induziert die einfachste solche Metrik verschwin- 
dende Kurvenlangen-Distanz am Raum der Figuren. Diese Entdeckung von Mi- 
chor und Mumford war der Ausgangspunkt fiir die Suche nach starkeren Metri- 
ken, die es ermoglichen sollten, Figuren anhand von fur die jeweilige Anwendung 
wichtigen Merkmalen zu unterscheiden. Sobolev-Metriken sind ein vielverspre- 
chender Ansatz dazu. Es gibt sie in zwei Varianten: Sogenannte auBere Me- 
triken, die von Metriken auf der DifFeomorphismcn-Gruppe des umgcbenden 
Raums induziert werden, und innere Metriken, die intrinsisch zur Figur defi- 
niert sind. In dieser Arbeit werden innere Sobolev-Metriken in einem sehr all- 
gemeinen Rahmen cntwickelt und prasc!ntic;rt. Es gibt kcinc Einscliriinkung auf 
die Dimension des eingebetteten oder umgebenden Raums, und der umgebende 
Raum muss nicht flach sein. Es wird gezeigt, dass die Kurvenlangen-Distanz 
von inncrcn Sobolev Metriken am Raum der Figuren nicht vcrschwindct. Die 
Geodatengleichung und die durch Symmetrien erzeugten Erhaltungsgr5fien wer- 
den hergeleitet, und die Wohldefiniertheit der Geodatengleichung wird gezeigt. 
Beispiele von numerischen Losungen der Geodatengleichung schliefien die Arbeit 
ab. 
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